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1.  INTRODUCTION 

The  smoke  chamber,  which  is  situated  at  MRL,  was  built  and  developed  to 
investigate  the  transmission  properties  of  current  and  novel  obscurants.  The 
chamber  can  be  used  to  investigate  transmission  in  the  visible  and  infrared, 
up  to  15  pm. 

This  manual  should  be  used  as  the  current  working  and  maintenance  manual;  it 
contains  sufficient  detail  to  duplicate,  set  up  and  maintain  the  facility. 

Earlier  reports  describe  the  associated  data  logging  system(ref . 1)  and  the 
philosophy  behind  the  smoke  chamber (ref . 2) . 


2.  HISTORY 

This  transmissometer  system  was  developed  in  1979  in  ERL  (IOC  Group)  as  its 
contribution  to  the  Australian  smoke  programme.  MRL's  contribution  is  the 
running  of  the  Smoke  Chamber  and  testing  and  development  of  various 
obscurants . 

The  visible  and  short  waveband  IR  region  was  measured  by  using  several 
discrete  passband  filters  switched  into  the  optical  train  in  turn,  and  the  IR 
band  from  2  to  15  pm  by  using  a  continuously  rotated  Circular  Variable  Filter 
(CVF) . 

For  the  IR  system,  since  low  level  signals  were  contemplated,  synchronous 
detection  was  chosen  as  desirable,  and  because  of  the  wide  spectrum 
(2  to  15  pm)  a  Golay  cell  was  required  as  the  detector,  this  being  the  only 
available  detector  to  conform  to  low  level  signals,  wide  dynamic  range  and 
wide  spectral  response. 

The  IR  source  needed  to  have  a  large  proportion  of  IR  content,  and  so  the 
choice  was  between  a  glow  bar  from  Oriel,  or  a  Nernst  Glower  (NG) .  A  NG  was 
chosen  because  of  its  higher  temperature  output  above  1000°K.  Two  suppliers 
of  NGs  are  Perkin  Elmer  and  Barr  and  Stroud  (used  in  the  Beckman  IR4 
Spectrometer) . 

At  the  time  of  system  design,  an  Oriel  Circular  Variable  Filter  set  (CVF)  was 
available  with  the  required  accuracy  in  wavelength  and  transmission  and  was 
used  to  design  a  three  filter  set  covering  a  range  of  2  to  15  pm  in  the  three 
sections : 

(a)  2.46  to  4.48  pm; 

(b)  4.47  to  8.04  pm;  and 

(c)  8.24  to  14.57  pm. 

Initially  this  system  was  designed  to  use  a  chart  recorder  and  the  appropriate 
data  reduction  carried  out  at  a  selected  number  of  points,  with  the  idea  that 
in  the  future  a  computer  data  logging  system  would  be  used:  hence  the  CVF  has 
a  graticule  marked  in  100  pulses/marks/segment,  but  only  every  one-tenth  could 
be  used  by  the  chart  recorder. 

The  concept  was  to  note  'no  input'  reading  level  from  the  Golay  cell,  with  the 
NG  blanked  off.  The  system  and  chart  recorder  are  then  started  and  the 
reading  level  recorded  with  no  smoke  present  ('no  smoke'  run).  Smoke  or  other 
obscurant  is  then  introduced  into  the  chamber  and  a  set  of  measurements  called 
the  'smoke'  run  are  recorded.  As  noise  is  a  problem,  several  sets  of 
readings,  ie  several  revolutions  of  the  CVF,  need  to  be  taken  and  each  point 
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(ie  each  wavelength  interval)  averaged  for  both  'no  smoke'  and  'smoke'  runs. 
There  are  usually  3  and  sometimes  5  revolutions  of  the  CVF  to  be  averaged  in 
this  way. 

2.1  The  electronic  design 

A  standard  laboratory  synchronous  detector  circuit  was  used,  taking  its 
reference  from  the  Oriel  chopper  at  the  Nernst  Glower  end  of  the 
transmissometer ,  and  generating  an  ac  signal  at  approximately  15  Hz.  A 
15  Hz  chopping  frequency  was  chosen  since  the  Go lay  cell  cannot  adequately 
follow  signals  faster  than  this,  and  still  retain  sufficient  gain.  This 
factor  dictated  the  speed  at  which  the  CVF  could  revolve. 

The  system  wavelength  resolution  is  defined  by  the  size  of  the  image 
focussed  on  the  CVF,  the  limited  speed  of  response  of  the  Golay  cell,  and 
the  scan  rate  of  the  filter.  Sampling  at  the  rate  of  458  samples/complete 
revolution  in  33.3  s  provided  a  resolution  of  100  samples/filter  segment, 
in  a  time  sufficiently  short  to  prevent  significant  changes  in  the 
obscurant  sample,  such  as  settling(ref . 1) . 

A  mass  produced  IR  source  and  detector  were  available  from  Hewlett  Packard 
as  an  integral  unit,  with  quite  high  resolution  and  speed.  This  HEDS1000 
was  used  for  reading  the  CVF  graticule  using  a  comparator  IC  LM311  with  a 
small  amount  of  hysteresis.  The  signal  was,  however,  too  fast  for  the 
chart  recorder  to  follow  and  so  each  100  pulse  train  was  divided  by  10, 
using  a  time-out/reset  circuit  which  incorporated  a  gated  free  running 
oscillator  into  a  counter,  the  output  of  which  triggered  the 
'synchronisation  pulse  detected'  monostable  which  is  set  for  1.5  s  pulse. 
The  same  triggering  action  resets  the  divide  by  10  counter. 

The  signal  from  the  Golay  cell  was  ac  amplified  through  a  low  noise 
amplifier  with  a  gain  of  5,  before  being  synchronously  detected;  some  gain 
and  isolation  being  necessary  to  get  sufficient  signal  to  be  handled  by  the 
synchronous  detector,  and  prevent  loading  of  the  high  impedance  Golay 
output  signal  by  periodic  clamping  of  the  synchronous  detector, 
respectively.  This  dc  signal  from  the  integrating  synchronous  detector 
was  passed  down  a  transmission  line  to  the  chart  recorder.  Since  the 
chart  recorder  could  not  follow  signals  greater  than  a  couple  of  cycles, 
this  system  appeared  quite  satisfactory. 

Deflections  corresponding  to  transmission  at  each  wavelength  of  concern 
were  measured  from  the  chart  records  and  tabulated  for  3  runs,  and  then 
averaged.  A  'no  smoke'  or  clear  air  run,  and  a  'smoke'  or  obscurant  run 
were  manipulated  in  this  way  before  the  mass  extinction  coefficient  for 
each  wavelength  could  be  calculated.  Since  the  mass  extinction  coefficient 
is  calculated  from  ratios  of  'smoke'  and  ’no  smoke'  readings,  absolute 
readings  were  unimportant,  as  long  as  the  signal  was  greater  than  the  noise 
and  precautions  were  taken  to  prevent  amplifier  'bottoming'  during  current 
alignment.  Therefore  long  term  deterioration  of  the  IR  source  or  the 
detection  system  was  basically  unimportant:  short  term  drift,  however, 
could  cause  errors. 

To  minimise  these  errors  the  Nernst  Glower  was  set  to  run  at  a  (low)  colour 
temperature  (1000°K),  and  for  stability  a  dc  power  supply  was  designed 
using  the  constant  current  collector  characteristics  of  the  transistor. 
Since  the  original  Nernst  Glowers  were  sintered,  switching  on  and  off 
resulted  in  a  high  failure  rate,  and  so  when  not  in  use  the  Nernst  Glowers 
were  put  in  ' standby 'mode,  using  sufficient  current  to  keep  them  self 
sustaining.  This  meant  a  running  current  of  0.7  A  and  a  standby  current  of 
0.2  A  for  the  3  mm  diameter  Nernst  Glower. 
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The  original  Nernst  Glowers  used  were  probably  spares  from  the  Beckman  IR4 
spectrometer  of  1.5  mm  diameter:  these  required  0.6  A  to  run  at  2100°K, 
and  0.2  A  for  'standby'.  These  Nernst  Glowers  were  superseded  by  the 
Perkin  Elmer  type  No  221-0451,  a  true  ceramic,  which  is  30  mm  long  and  has 
a  diameter  of  3  mm,  requiring  similar  currents  able  to  run  at  1000°K. 

Both  Nernst  Glowers  require  heating  to  start  them,  and  until  recently 
heaters  from  the  Beckman  IR4  were  used  -  these  required  some  5  A  at  10  V. 

The  Nernst  Glower  housing  and  lens  assembly  were  modified  Oriel  types 
Nos  6361  and  6362  respectively.  Mica  brackets  held  the  'glower  and 
heaters  in  position,  and  not  being  mechanically  reliable  needed  mica  wedges 
to  separate  the  elements  to  their  desired  positions  subsequent  to  any 
movement  and/or  installation. 

A  Nernst  Glower  is  a  ceramic  mixture  of  zirconium,  yttrium,  thorium  and 
other  oxides,  which  is  heated  externally  to  400°C  before  it  will  pass 
sufficient  current  to  sustain  itself. 

The  Smoke  Chamber  system  w^as  computerised  using  an  LSI  11/03  computer;  the 
PDP  11/34  mainframe  was  used  to  develop  the  programmes  which  were  then  run 
on  the  LSI  11/03.  A  configuration  was  devised  to  conform  to  the  initial 
specifications,  and  was  initially  oriented  to  the  2  to  15  pm  region  which 
is  still  the  region  of  main  area  of  interest.  After  installation,  a 
'setting  up'  procedure  needed  to  be  closely  followed.  The  initial 
extinction  coefficient  results  were  no  better  than  those  obtained  using  the 
10  points  measured  from  the  chart  recorder,  but  subsequently  the  system  was 
improved  to  such  an  extent  that  the  test  results  from  various  calibration 
filters  eventually  conformed  remarkably  well  with  published  data,  as 
various  problems  were  identified  and  solved. 

2.2  Future  developments 

Since  the  electronics  are  spread  over  several  boxes  and  parts  of  the  system 
are  duplicated  and  or  redundant  a  tidier  electronic  design,  which  is 
functionally  identical  to  the  existing  system,  is  contemplated.  In  the 
event  of  a  catastrophic  failure,  it  is  desirable  that  the  existing  system 
can  be  readily  replaced  by  this  updated  version. 

Also,  since  the  Golay  cell  is  expensive  and  has  a  comparatively  slow 
response  (25  Hz),  a  faster  detector  system  should  be  found,  which  would 
then  increase  the  speed  of  data  logging  by  at  least  ten  times.  This  should 
make  the  system  perform  closer  to  a  real  time  situation  ie  less  probability 
of  particle  differential  fall-out  due  to  the  33  s  CVF  revolution. 

Another  point  to  be  considered  is  the  heating  of  the  Platimum  NG  heaters 
which  will  probably  require  a  constant  current  source  because  of  the  wide 
tolerances  of  the  CZ12  NTC  resistors  used:  ie  their  heating  time  constants 
appear  to  change  by  at  least  3:1,  which  causes  problems. 


3.  PRESENT  SYSTEM  DESCRIPTION 


3 . 1  Source 

The  Nernst  Glower  housing  is  still  being  developed.  The  original  housing 
was  one  from  the  Oriel  range  but  modified  to  hold  the  Nernst  Glower.  A  new 
chopper  which  was  installed  in  the  side  wall  caused  noise  and  vibration  so 
this  was  changed  to  a  vertical  mounting.  (Vibration  will  reduce  the  life 
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of  the  Nernst  Glower  and  so  resonances  had  to  be  eliminated  by  muffling  or 
design.)  A  5  mm  square  aperture  is  placed  at  the  focal  point  of  a  bloomed 
germanium  transmitting  lens.  (See  figure  1). 

A  Nernst  Glower  is  designed  to  operate  in  air  at  a  temperature  of 
2100  K  (1827°C).  What  appears  to  be  a  polarising  problem  within  the 
Nernst  Glower,  when  using  a  dc  power  source  beyond  0.8  A,  is  eliminated  by 
using  a  controlled  ac  power  supply.  This  ac  power  supply  was 
designed  to  drive  the  Nernst  Glower  from  its  "dc  @  1000°K"  capability  to  an 
upgraded  limit  of  2000°K.  This  alone  should  give  a  4:1  increase  in  useful 
IR  energy  in  the  range  2  to  15  pm. 

The  Nernst  Glower  has  a  negative  resistance  coefficient  with  a  slope 
-7.3  ft/A  at  an  operating  temperature  of  2100°K,  the  current  required  being 
2.1  A.  Therefore  the  Nernst  Glower  power  supply  needs  to  deliver  92  W, 
but  to  first  strike  the  Nernst  Glower  a  voltage  in  excess  of  70  V  is 
required.  (See  figure  2). 

A  Cadmium  Sulphide  (CdS)  detector  is  used  as  the  sensing  element  for 
luminous  flux  control.  This  detector  is  in  a  bridge  circuit  with  the 
^NG  potentiometer  feeding  into  a  high  gain  operational  amplifier 

(Ay  =  100),  then  to  the  controlling  transistor  network  to  the  magnetic 

amplifier.  Circuit  layout  problems  have  been  taken  into  consideration  (ie 
earth  loops  considered) . 

The  power  supply  is  a  controlled  current  unit  using  a  magnetic  amplifier, 
which  was  designed  in  AEL,  for  control  from  a  single  transistor.  This 
magnetic  amplifier  has  been  designed  to  produce  a  near  sine-wave  output 
even  when  operated  at  maximum  dc  input  (which  provides  maximum  power  to  the 
NG) .  The  objective  of  the  smoothly  varying  ac  waveform  was  to  minimise 
rapid  current  changes  and  to  prevent  build  up  of  polarising  effects  in  the 
NG. 


Another  problem  which  is  not  normally  encountered  is  the  heaters  required 
to  start  the  NG.  Since  the  NG  runs  at  such  a  high  temperature,  platinum 
heaters  were  the  choice  of  the  NG  manufacturers .  This  material  has  a  large 
positive  temperature  coefficient  such  that  at  the  running  temperature  it 
requires  60  V  at  2  A  (30  Q)  but  when  cold  the  resistance  is  much  lower  and 
it  draws  excessive  current  from  the  transformer.  Therefore  a  Negative 
Temperature  Coefficient  (NTC)  resistor  had  to  be  found  that  could  handle 
2  A  while  running;  its  resistance  has  to  drop  to  parts  of  ohms  when  hot, 
and  be  at  least  30  Q  resistance  when  cold.  A  CZ12  was  the  only  device 
available. 

Once  the  NG  became  self  sustaining  the  heaters  needed  to  be  switched  off, 

and  so  a  unique  method  was  devised  to  achieve  this.  Using  the  forward 

voltage  drops  across  3  diodes  (1N5626),  and  a  small  bridge  rectifier  and 

smoothing  capacitor,  to  give  a  dc  'signal';  an  opto  isolator  of  the  type 

H13A1  was  used  to  sense  the  NG  current.  The  sensed  I  is  fed  into  a 

NG 

comparator  which  is  set  to  switch  at  I  greater  than  0.3  A,  which  is  the 

minimum  sustaining  current  for  the  NG.  A  solid  state  relay  is  used  to 
switch  off  the  heaters  -  this  method  cannot  burn  contacts  and  is  a  zero 
voltage  crossing  switch  which  reduces  RFI  potential  problems. 

3 . 2  Chopper 

A  controlled  chopper  was  decided  to  be  necessary  and  a  loudspeaker  was 
found  to  have  the  appropriate  characteristics.  A  12  W  8  ft  high  resilience 
loudspeaker  has  the  necessary  depth  of  throw  and  frequency  response  with 
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the  power  to  drive  a  chopper  blade  over  a  7  mm  distance.  The  shaft  should 
be  longitudinally  polished  to  reduce  wearing  of  the  bearings,  and  the 
chopper  positioned  vertically,  also  to  reduce  bearing  wear,  but  more 
importantly  to  reduce  vibration  and  noise.  (See  figures  1  and  3). 

The  chopper  driver  consists  of  two  transistors,  one  of  which  biases  the 
chopper  blade  from  half  open  to  fully  open,  and  the  other  drives  the  blade 
fully  closed,  on  command. 

Adjustments  for  various  choppers  (resilience  change)  is  controlled  by 
either  parallelling  27  Q.  12  W  resistors  to  adjust  the  'opening'  bias, 
and/or  to  series  1  2  12  W  resistors  for  driving  the  blade  closed.  Normally 
this  circuit  and  chopper  combination  will  function  up  to  50  Hz,  and  some  as 
high  as  70  Hz:  the  requirement,  in  fact,  is  only  15  Hz  (13.7  Hz).  No 
other  device  so  far  has  been  found  to  give  at  least  5  mm  throw  and  at  least 
20  Hz  response. 

3.3  Sampling  and  detection 

The  radiometer  consists  of  various  parts  -  some  mechanical,  others  optical 
and  electronics  (see  figures  A  to  7,  Table  1). 

The  electronics  block  diagram  (and  schedule)  are  shown  in  figure  8  and 
Table  2. 

IR  radiation  from  the  NG  crosses  a  5  m  path  length  and  enters  a  germanium 
lensed  telescope  focusing  on  a  slit  1.2  mm  wide  with  the  CVF  wheel  directly 
behind.  The  resolution  varies  with  the  filter  segment  in  use  and  is 

1.99  x  10~2  ym,  3.67  x  10  2  ym,  6.4  x  10  2  ym  respectively,  as  resolved  by 
the  computer (ref . 1,2) . 

Being  a  pneumatic  operating  device,  and  therefore  its  microphonic  nature, 
the  Golay  cell  is  mounted  on  an  antivibration  mount.  It  is  positioned 
behind  the  CVF  (see  figures  9  and  10). 

The  chopped  IR  energy,  being  generated  by  the  'driver  chopper',  is 
amplified  within  the  detector  housing  before  being  cabled  to  another  low 
noise  ac  coupled  amplifier  which  feeds  into  the  synchronous  detector 
system. 

The  CVF  is  driven  via  a  gearbox  from  a  mains  synchronous  motor  and  the  CVF 
wheel  takes  33.3  s  to  turn  one  revolution.  The  CVF  is  a  3  segmented 
filter  set,  set  in  a  wheel  to  which  is  attached  a  glass  rim  with  a  vacuum 
deposited  nichrome  graticule  which  has  100  marks / s egment ,  and  a 
synchronising  mark  just  prior  to  segment  No  1  (see  figures  11  and  12). 
These  marks  are  detected  using  a  HEDS1000  light  and  detector  unit  (see 
figures  13  and  14)  which  has  a  sharp  focus  to  detect  each  mark  as  it  passes 
through  its  operating  field.  The  signal  is  'squared  up'  using  an  LM311 
comparator,  the  output  of  which  is  termed  'Read  Pulses'.  (The  HEDS1000 
was  originally  designed  by  Hewlett  Packard  to  be  used  in  bar  code  readers). 

These  'Read  Pulses'  go  to  four  areas,  one  being  the  synchronous  detector; 
another  the  data  logger,  the  third  area  generates  a  more  suitable  pulse 
train  for  writing  on  the  backup  chart  recorder,  and  the  fourth  area 
generates  suitable  delays  and  pulse  widths  for  the  controlled  chopper  (and 
the  synchronous  detector)  (see  figure  16) .  The  chart  recorder  can  only 
resolve  one-tenth  the  'Read  Pulses'  speed  and  so  a  divide  by  ten  circuit  is 
used,  with  a  pulse  stretching  monostable  generating  a  pulse  suitable  for 
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use  with  the  chart  recorder.  This  circuit  has  a  gated  (by  the 
Read  Pulses)  oscillator  and  counter  which  gives  a  'carry  out'  signal  only 
on  the  synchronising  pulse  (see  figures  15  and  26). 

The  'Read  Pulses'  from  the  HEDS1000  comparator  latch  the  synchronously 
detected  IR  signal  into  a  sample  and  hold  circuit  (see  figure  16).  This 
had  been  found  necessary  since  the  synchronised  signal  is  only  valid  at  one 
point  in  time,  due  to  the  Golay  cell  response  and  the  synchronous  detector 
type  used,  so  the  chart  recorder  printout  is  basically  a  stepped  trace. 

It  should  be  noted  that  electrically,  the  Read  Pulses  are  shifted  one  pulse 
which  would  cause  the  apparent  CVF  position  marker  to  be  out  of  step  with 
the  IR  signal;  but  since  only  96  of  the  100  pulses  available  are  used,  the 
software  takes  this  apparent  discrepancy  into  account  by  only  using,  and 
aligning  the  correct  wavelength  position  of  the  CVF.  Hence  the  software 
looks  at  100  pulses,  but  ignores  one  pulse  at  each  end  because  of  non 
linearities  at  the  filters  ends.  The  first  pulse,  however,  is  used  to 
initiate  the  chopping  sequence.  The  chart  recorder  cannot  resolve  less 
than  4  pulse  widths  and  so  was  not  affected  by  this  problem  in  the  past. 

Timing  and  phasing  of  the  various  driver  signals  derived  from  the 
CVF  Read  Pulses  need  to  be  tuned  for  maximum  signal  from  the  Golay  cell, 
and  so  are  probably  more  complex  than  need  be. 

The  timing  is  initiated  from  the  CVF  Read  Pulses  HEDS1000  comparator 
circuit.  As  it  has  been  said  before,  it  branches  to  four  main  areas:  the 
one  considered  now  is  to  the  synchronous  detector  circuitry  (see  figures  17 
to  21).  This  has  two  delay  branches,  one  delay  approximately  30  ms 
triggers  the  chopper  driver  at  the  NG  end;  the  other  delays  and  'latches' 
the  synchronous  detector  -  an  initial  delay  of  approximately  5.5  ms  and 
then  a  31  ms  'on'  time  for  the  synchronous  detector. 

These  three  delays  (monostables)  are  tuned  so  the  computer  does  an  analogue 
to  digital  conversion  at  the  peak  of  the  IR  signal,  (see  figure  22).  The 
top  is  comparatively  flat,  and  the  edges  can  and  do  vary  in  time  due  to  the 
variation  in  spacing  of  the  CVF  graticule  marks. 

3.4  Power  supplies 

The  power  supplies  need  not  be  described  in  any  detail  as  they  are  fairly 
standard,  mostly  of  the  IC  type.  However,  the  Golay  cell  power  supply 
(see  figures  23  to  25)  can  be  classified  as  a  high  loss  (high  impedance) 
type  which  relies  on  a  high  voltage  (approximately  130  V  dc)  and  large 
resistors  (18  kft)  for  filtering,  and  the  output  only  becomes  24  V  when 
loaded  with  the  appropriate  Golay  amplifier  circuit.  That  is,  the 
filtering  in.  particular  divides  down  the  ripple  content.  The  Golay' s 
'lamp'  power  supply,  however,  is  a  discrete  component  circuit.  This  lamp 
circuit  also  incorporates  a  constant  current  device  (within  the  Golay  cell 
housing) . 


4.  COMPUTER  INTERFACING  AND  CONFIGURATION 
4.1  Computer  system  hardware 

The  computer  data  logger  is  an  LSI  11/03  Q-bus  frame  and  CPU,  with  a 
parallel  card,  two  serial  cards,  and  analogue  to  digital  conversion  card 
DRV 1 1 ,  DLV11 ,  DLV11-J,  DT2764  respectively,  talking  to  the  outside  world 
(see  figure  27).  The  serial  cards  DLV11  and  DLV11-J  talk  to  the  console 
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(keyboard  and  VDU) ,  plotter  and  printer,  all  through  software  'handlers' 
and  so  are  not  normally  associated  directly  with  user  control.  DEC 
standard  addressing  is  normally  used. 

Because  'bit  15'  can  be  easily  read  by  the  Macro  instruction  TST  the 
parallel  card  'DRV11  bit  15  input'  is  the  logical  bit  to  use  to  fast  read 
the  pulses  from  the  CVF  Read  Pulses. 

The  DT2764  is  a  12  bit  analogue  to  digital  converter  card  specifically 
designed  by  Data  Translation  to  suit  the  DEC/LSI  Q-bus .  It  has  been 
configured  for  single  ended  input,  which  then  dictates  that  16  channels  are 
available  (the  first  channel  is  used),  and  an  input  voltage  range  of  ±10  V, 

using  2S  compliment  at  the  output.  This  gives  the  ADC  output  to  input 
range  of  2047  (decimal)  equivalent  to  9.99  V,  and  -2048  (decimal) 
equivalent  to  -10.0  V. 

It  should  be  noted  that  the  computer  must  be  switched  on  before  the  Golay 
cell  radiometer  output  becomes  meaningful,  as  the  computer,  when  switched 
off,  loads  the  CVF  Read  Pulses,  and  prevents  them  from  being  read  even  by 
the  chart  recorder. 

4.2  Chart  recorder  backup 

The  chart  recorder  used  is  a  Curken  model  250-3A.  The  radiometer 
electronics  produces  the  'divide  by  10  Read  Pulses'  which  use  channel  No  1 
of  this  3-pen  X-T  recorder  to  display  the  'position'  of  the  IR  on  the  CVF. 
A  trace  of  about  20  mm  high  is  sufficient,  and  this  allows  the  majority  of 
paper  width  to  be  used  for  both  IR  and  visible  radiometer  results. 
Channel  No  3  records  the  visible  radiometer  which  has  ten  filter  positions, 
but  three  are  used,  namely,  eye-response,  1.06  pm,  and  0.72  ym  filters; 
while  channel  No  2  shows  the  IR  vs  wavelength  response.  The  visible 
radiometer  system  has  not  yet  been  incorporated  in  the  computer  system. 

These  three  signals  should  be  carried  from  the  radiometer  table  at  the  end 
of  the  smoke  chamber  to  the  computer  room,  which  also  houses  this  chart 
recorder,  via  twin  shielded  cables,  to  reduce  interference  both  into  and 
out  of  the  cables.  Earth  loop  problems  are  reduced  by  referring  the 
signals  back  to  the  radiometers,  and  only  earthing  the  shield  at  one  place. 
A  10  kfi  or  higher  resistor  could  be  taken  from  signal  return  to  earth  to 
prevent  signals  floating. 


5.  SETTING  UP/TUNING  THE  SYSTEM 

To  obtain  the  best  possible  performance  from  this  overall  system,  each  stage 
needs  to  be  individually  optimised  and  then  collectively  tuned.  Using  the 
proposed  methods,  an  increase  of  perhaps  five  times  in  sensitivity  could  be 
achieved,  compared  to  only  trimming  the  sections  individually.  The 
radiometer  only  'takes'  relative  readings  and  so  overall  trimming  improves  the 
noise  figure  and  dynamic  range;  the  accuracy  is  basically  unaffected. 

5.1  Nernst  Glower 

The  Nernst  Glower  needs  to  be  run  as  hot  as  possible  without  melting  its 
connecting  leads.  The  top  lead  is  subjected  to  more  heat  stress  via 
radiation  and  convection  than  the  bottom  one,  and  so  an  empirical  method  is 
used  to  determine  maximum  Nernst  glower  temperature  (see  figure  28). 
Platimum,  the  lead  material,  melts  at  1942°K  but  is  approximately  3  mm  away 
from  the  Nernst  glower  which  could  therefore  be  run  200°K  higher.  The 
better  scheme  would  be  to  keep  the  leads  well  below  their  melting  point  and 
have  the  Nernst  Glower  temperature  set  to  about  2000°K  or  1730°C  (which  can 
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be  read  on  an  optical  pyrometer).  For  this  temperature  to  be  reached  the 
Nernst  Glower  needs  a  current  of  1.8  A.  With  efficient  reflectors  ’this 
current  may  be  reduced,  and  possibly  a  decrease  in  temperature  could  be 
tolerated.  A  four  times  improvement  in  IR  energy  is  sought  compared  with 
the  original  dc  drive  and  1000°K  response. 

To  adjust  the  Nernst  Glower  power  supply  (see  figures  29  to  32):  set  the 
run/standby  switch  to  'run',  set  the  'set  standby  current'  potentiometer 
to  maximum  resistance,  set  the  'heater  off’  potentiometer  to  the  'high* 
side  (positive  end),  and  set  the  Nernst  Glower  current  potentiometer  to 
approximately  1  kfl  (one-fifth  of  maximum  resistance).  Place  a  voltmeter 
across  the  heater  transformer  primary,  setting  the  meter  for  250  V  ac 
range.  On  switching  the  mains,  the  'mains'  indicator  should  glow,  and  no 
volts  read  on  the  meter.  Turn  the  'heaters  off'  potentiometer  down  until 
the  meter  reads  mains  voltage  of  230  V  ac.  The  heaters  will  take  about 
20  s  to  begin  to  glow,  and  heat  the  Nernst  Glower.  This  delay  is  due  to 
the  CZ12  negative  temperature  coefficient  resistance  (NTC)  which  prevents 
current  surge  into  the  platinum  heaters.  After  about  90  s  the 

Nernst  Glower  comes  up  to  temperature  as  set  by  the  I  potentiometer: 

• 

there  is  normally  an  overshoot  of  perhaps  0.3  A  for  a  second  before 
dropping  back  to  its  set  running  current.  The  heaters  should  also  have 
switched  off  (ac  voltmeter  reading) .  The  Nernst  Glower  current  should  be 
about  1.4  A.  ^  Switch  to  'standby',  and  adjust  standby  current  to  0.3  A  with 
the  standby  potentiometer:  the  heaters  should  not  switch  on  at  this 
level.  Reduce  the  'set  I^g '  potentiometer  so  ING  drops  to  0.2  A  and  then 

adjust  the  'set  heaters  off'  potentiometer  until  the  heaters  switch  on  (ac 
volt  meter).  Before  the  heaters  begin  to  glow,  readjust  the  I 

potentiometer  back  to  0.3  A,  whereby  the  heaters  should  again  switch  off. 
Put  the  power  supply  into  the  'RUN'  mode  and  readjust  the  INQ  potentiometer 

to  set  the  current  to  1.6  A.  (This  reduces  the  IR  level  somewhat  but 
increases  the  Nernst  Glower's  life  -  unless  more  IR  is  really  required). 

This  power  supply  should  then  be  ready  to  use. 

The  Nernst  Glower  bar  or  heater  wires  should  never  be  touched  because 
contamination  can  destroy  these  parts. 

5 . 2  Chopper 

The  chopper  needs  to  be  checked  for  proper  operation.  With  the 
Nernst  Glower  lens  assembly  removed  (the  85  x  60  mm  flanged  tube),  look  at 
the  chopper  blade  to  determine  its  rest  position  (see  figure  1).  It  should 
be  approximately  half  way  across  the  5  mm  square  aperture.  Switch  on  the 
NG  power  supply  (with  chopper  connected)  and  note  that  the  chopper  blade 
should  pull  down  and  fully  open  the  aperture.  Drive  the  chopper  input 
terminals  with  a  0  to  +5  V  square  wave  at  15  Hz,  and  note  the  swing  of  the 
chopper  blade.  It  should  fully  close  off  the  aperture.  There  are  two 
points  to  watch  while  trimming  this: 

(1)  the  blade  should  have  an  overall  travel  of  6  mm  minimum  and 
completely  open  and  completely  close  the  aperture;  and 

(2)  the  chopper  should  not  bottom  on  its  housing  in  either  direction. 

If  the  chopper  travel  is  critical  then  both  mechanical  and  electrical 
adjustments  are  required.  Mechanical  adjustments  could  take  the  form  of 
shims  under  the  chopper  housing,  or  replacing  the  blades  (and  this  is  not 
recommended).  Electrical  adjustments  are  simpler. 
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Without  any  square  wave  signal  the  chopper  should  open  the  aperture; 
however,  if  it  closes  then  the  chopper  drive  leads  are  crossed.  If  the 
opening  travel  is  not  far  enough  (or  too  far)  adjust/pad  the  27  Q 
resistors.  When  driving  the  chopper  circuit  with  a  15  Hz  5  V  signal 
adjust  the  1  Q,  resistors  to  effect  sufficient  travel  of  the  blade. 

5.3  Source  collimator  lens 

The  lens  assembly  is  replaced,  and  the  NG  and  lens  pointed  such  that  at 
least  a  clear  3  m  path  length  is  needed  for  collimating  the  IR  beam.  The 
area  should  be  darkened  to  IR  signals  so  that  stray  signals  will  not  be 
read  by  the  aligning  IR  detector,  and  mark  the  receiving  area  to  show 
limiting  (and  centre)  positions  of  view  of  the  NG  assembly.  The  objective 
is  to  get  as  small  and  uniform  a  spot  as  possible  which  is  both  round  and 
has  the  IR  level  consistently  flat  over  that  area,  by  adjusting  the  focus 
of  the  lens.  A  hand  held  PbS  detector  circuit  is  suitable  as  an  indicator 
for  this  procedure.  Refer  figure  33.  The  IR  beam  axis  should  be 
determined  and  a  line  marked  on  the  source  housing  to  aid  alignment  with 
the  Golay  radiometer's  telescope.  This  position  is  trimmed  later  when 
overall  system  tuning  is  carried  out. 

5.4  Sampling  indicator  (read  pulses) 

If  the  CVF  HEDS1000  detector  needs  to  be  shifted/ adjusted,  then  the 
computer  must  be  recalibrated.  Focusing  of  the  HEDS1000  detector  needs 
care.  Measurements  are  taken  from  the  top  face  of  its  housing  and  are 
typically  4.27  mm  to  target  surface.  There  may  be  a  variation  of  0.5  mm 
for  maximum  signal  for  best  focus,  and  so  the  focus  distance  should  only 
have  a  ±0.2  mm  discrepancy  from  its  natural  focus;  its  image  diameter  is 
0.17  mm.  The  central  wavelength  of  the  detector  is  700  nm.  There  is  no 
adjustment  for  the  HEDS1000  after  it  is  soldered  into  position,  so  this 
information  must  be  implemented  precisely  for  aligning  (installing)  a  new 
detector.  In  service,  using  the  glass  substrate  with  nichrome  markings 
(see  figure  12),  its  output  to  the  LM311  comparator  should  be  in  the  order 
of  0.9  V  PP  with  an  average/dc  level  of  +1.4  V  dc  (see  figure  13).  This 

1.4  V  dc  is  also  the  level  at  which  the  comparator  is  set.  'No  signal'  is 
equivalent  to  'low'  V  out.  This  signal  is  available  at  the  rear  end  of 
the  Golay  cell  housing  on  a  BNC  connector  marked  'Read  Pulse'  (see 
figure  6). 

(a)  The  CVF  itself  houses  three  AX  segments  (see  figure  11);  however. 
Nos  2  and  3  sections  have  been  installed  back  to  front  and  only  give  the 
pretence  of  normality  when  the  results  are  plotted.  This  is  corrected 
in  the  computer  programme,  but  only  needs  to  be  recognised  as  a 
phenomenon  here . 

To  adjust  any  of  the  Golay' s  signal  path  the  CVF  needs  to  be  put  in  its 
maximum  transmission  position,  which  is  200°  as  marked  on  the  position 
graticule  as  viewed  through  the  port  above  the  IOC  Group's  Logo,  and  the 
drive  motor  to  the  CVF  switched  off.  The  'Read  Pulses'  lead  is  taken  off 
the  Radiometer,  from  the  synchronous  detector,  and  plugged  into  the 
'CVF  Simulator'  box  (see  figure  34)*.  This  then  gives  a  constant  output 
from  the  Golay  cell  that  for  all  intents  and  purposes  drives  the  chopper, 
synchronous  detector  and  computer  as  though  the  system  is  fully  functional. 
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5.5  Golay  detector 

Since  the  Golay  cell  box  should  never  be  opened,  only  a  superficial 
escription  needs  to  be  given.  The  Golay  cell  is  a  pneumatic  cell  that 
changes  its  volume  when  energy  falls  onto  its  input  window;  the  cell  wall 
on  the  aft  side  moves  and  deflects  a  light  beam  from  a  constant  current 
light  emitting  device.  There  are  two  sections  to  the  Golay  'electronics' 
one  being  the  light  emitter  and  detector,  and  the  other  is  the  low  noise  dc 
amplifier.  The  cell  itself,  of  the  Golay  cell,  comprises  two  sections 
with  a  slow  leak  between  them.  When  in  the  presence  of  a  high  ambient 
energy  source  the  working  part  of  the  cell  can  go  into  a  nonlinear  region, 
but  recovers,  so  the  built  in  preamplifier  will  also  work  over  its  linear 
range.  The  signal  of  concern  'rides'  on  the  ambient  light  level  The 
signals  to  date  are  in  the  order  of  6  mV  PP  of  random  noise  (no  IR  input! 
and  600  mV  PP  for  maximum  signal. 

5.6  Detector  telescope 

The  Golay  cell's  telescope  should  not  need  to  be  adjusted,  but  may  be 
cleaned  with  an  alcohol  dampened  soft  cloth,  being  careful  not  to  scratch 
the  germanium  lens  s  bloomed  surface.  If  need  be,  the  telescope  is 
focused  by  using  the  CVF  Simulator,  tuning  the  telescope  for  maximum  Golay 
cell  output.  Since  the  optical  and  mechanical  axes  are  not  coincident  if 
this  focus  is  varied,  then  the  whole  system  needs  to  be  recalibrated. 

To  ’point’  the  radiometer  for  maximum  signal  the  CVF  Simulator  needs  to  be 
used  while  the  CVF  itself  has  been  set  to  the  200°  mark  and  the  motor 
switched  off.  The  Golay  output  is  monitored  by  a  CRO,  and  the  whole 
housing  moved/wedged  etc  for  maximum  output  (the  objective  being  to  centre 
the  telescope  objective  lens  in  the  beam  from  the  Nernst  Glower,  and  to 
ensure  that  the  lens  axis  is  sufficiently  aligned  with  the  beam  axis  to 
maximise  signal).  The  electronics  can  be  tuned  for  maximum  output  from  the 
synchronous  detector. 

5.7  Electronic  circuit  tuning 

The  IR  signal  from  the  synchronous  detector  gets  a  final  dc  amplification 
with  an  appropriate  offset  voltage  to  force  the  analogue  signal  to  the 
computer  to  be  always  positive,  and  not  overloaded  by  varying  its  gain  (see 
figure  16).  Positive  9.8  V  should  be  maximum  signal  for  the  analogue  to 
digital  converter.  Set  the  offset  potentiometer  such  that  with  no  IR  input 
to  the  Golay  cell,  with  the  telescope  'capped',  the  output  to  the  A/D 
converter  should  now  be  small  but  positive  -  say  +15  mV.  This  should  be 
done  with  no  IR  signal,  but  the  CVF  itself  should  be  turning.  The  time  at 
which  the  computer  takes  a  reading  is  critical  (see  figure  22).  With  too 
short  a  time,  the  signal  is  still  rising  up  to  the  value,  but  with  too  long 
a  time  the  synchronous  detector  has  begun  to  'clamp',  preparing  for  the 
next  signal  pulse. 

In  order  to  generate  the  proper  timing  of  chopper  drive,  synchronous 
detection  and  sample  and  hold  circuits,  the  first  (of  the  100  pulses  in 
each  segment)  pulse  detected  initiates  the  timing  sequence  controlled  by 
monostable  delays  such  that  the  first  analogue  to  digital  conversion 
reading  is  invalid. 


This  area  of  the  circuit  operation  is  as  follows:  when  the  CVF  pulse  goes 
igh  a  Read  Pulse  is  generated  which  is  detected  by  the  computer  (and  the 
divide  by  ten  circuit  for  the  chart  recorder),  and  if  it  is  detected  as  a 
read  pulse  and  not  a  'synchronising  pulse'  an  A/D  conversion  is 
performed,  and  stored.  But  since  on  the  first  pulse  (of  each  segment)  the 
synchronous  detection  and  integration  circuitry  have  not  the  desired  time 
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scale,  this  reading  would  be  erroneous.  But  as  well,  since  this  first 
reading  also  appears  to  suffer  from  edge/end  effects  (reflections)  from  the 
mounting  (of  the  filter  segments)  and  out  of  tolerance  (linearity)  of  the 
wavelength  vs  angular  rotation  of  the  filter,  this  first  reading  is  set  to 
zero,  as  an  error,  and  therefore  not  graphed. 

The  circuit  is  such  that  a  reading  on  the  A/D  converter  is  taken  just  prior 
to  the  chopper  closing  off  the  IR  signal  at  the  Nernst  Glower 
(approximately  3  ms  chopper  reaction  time).  Timing  of  the  chopper  drive 
monostable  is  not  especially  critical,  but  is  a  compromise  of  maximum 
available  signal  being  effectively  read,  and  settling  time  for  the 
synchronous  detector  to  define  the  zero  point  of  the  signal.  Period  is 
approximately  73  ms,  with  the  chopper-shut  time  of  30  ms  (29  ms). 

The  synchronous  detector  needs  a  delay  from  this  reference  edge  of  4  ms  so 
the  IR  signal  can  be  read  by  the  signal  reading/ latching  system  before  the 
synchronous  detector  switches  over  to  admitting  the  IR  signal.  The 
synchronous  detector  passes  this  IR  signal  for  34.5  ms,  but  the  Golay  cell 
will  not  peak  until  about  the  26th  ms;  therefore  there  is  only  about  3  ms 
in  which  the  signal  out  of  the  last  amplifier  is  at  its  peak  and  is  there 
to  be  read  by  the  A/D  converter.  Therefore  the  'art'  of  tuning  this  system 
is  to  reduce  droop  at  point  a)  and  make  point  b)  as  high  and  flat  as 
possible  (on  the  sample  and  hold  (S/H)  input  waveform)  using  the  Read 
Pulses  on  the  reference  edge.  (Rising  edge  latches  the  S/H  circuit)  (see 
figure  22). 

5.8  Chart  recorder  time  reference 

Because  the  chart  recorder  can  only  follow  one-tenth  the  speed  of  the 
computer,  the  100  Read  Pulses  need  to  be  divided  down  to  10  modified 
Read  Pulses,  and  the  synchronising  pulse  needs  to  be  broadened.  The 
Read  Pulses  are  put  through  a  divide  by  10  counter  -  a  4017  -  which  then 
triggers  a  1  s  monostable.  The  'high'  value  of  the  Read  Pulses  enables 
the  340  Hz  oscillator  which  is  counted  by  the  4024.  Normally  the  count  is 
cut  short  by  the  falling  edge  of  the  Read  Pulse  and  resets  this  counter, 
but  when  the  synchronising  pulse  is  encountered,  the  4024  terminally 
counts;  it  resets  the  4017  and  triggers  a  1.5  s  monostable.  These  two 
monostable  pulses  are  summed  and  fed  to  the  chart  recorder  channel  No  1. 


6 .  PROGRAMMING 

Two  programmes  MRLV20  and  MRLGRP  form  a  suite  to  obtain  and  display  IR  data 
from  the  Smoke  Chamber . 

The  philosophy  is  to  keep  the  current  version  of  the  running  programme,  both 
MRLV20  and  MRLGRP  along  with  IRUN(=1) ,  on  the  'data  disc'  (DK:)  which  then 
ensures  that  the  data  gathered  is  relevant,  and  if  an  outdated  programme  is 
used  it  should  be  readily  detected  by  reading  the  onboard  (disc)  programme. 

The  programme,  MRLV20,  for  the  Data  Logger  section  was  mainly  written  in 
Fortran  for  ease  of  programme  maintenance,  and  is  Menu  driven  for  ease  of  use. 
Each  selection  is  a  subroutine  with  a  Common  Block  (of  memory)  for  passing 
data  and  parameters.  Reference  1  describes  the  data  reading  and  storing 
Macro  subroutine,  READ43.MAC,  with  sufficient  detail  to  not  need  describing 
here.  (The  listing  is  in  Appendix  I.) 

MRLGRP  uses  existing  data  on  disc,  and  was  written  to  graph  not  only  Mass 
Extinction  Coefficient  but  also  the  Raw  Data  ('smoke'  and  'no  smoke'  readings 
on  the  same  graph).  Also  it  provides  a  Percentage  graph  of  'smoke  run'  for 
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f±lter  calibration  results  for  checking  system  integrity  fa 
rfsuJt?  57  “  ^  3  Pr°blem  ^  be  S°lved  by  usin§  a  'lower  level* 

Umitations  (PDP  n/34,  VT103/440)  negate  a  single  programme 
editaMe  at  one  session,  these  two  programmes  needed  to  be  suitably  divided 
into  several  sections  and  linked'  (the  PDP  word  for  joining  and  completing 
compiling  of  several  programmes  to  make  one  large  one).  P  S 

6.1  Graphing  routines 

(D  .MRLGRp.F°R  (figures  35  to  38  show  flow  charts,  and  Appendix  II  the 
is ting) ,  is  the  filename  which  comprises  the  MAIN  routine  (figure  35) 
consisting  of  initialising  and  setting  of  data  blocks,  variables  and 
manipulating  the  Menu.  There  are  four  subroutine  call’s,  selected  iTa 

RAW  (IFTn°nh-  r  RERUN  *hich  reads  run  data  from  the  disc,  and  three  to 
^fL(pFIL)  Whl^h  PaSS  thS  aPPr°Priate  flag  for  the  required  graph  type 

End  Programme  waits  a  short  time  before  the  computer  returns  to^he 

Operating  System  and  presents  the  period  (.)  prompt. 

cnLentTlTv  “°LJN  .ri“atf  Une'>:  <^-36,  draws  to  the 
coordinates  (calculated  in  the  subroutine  RAW  (IFIL)) 

return.  ^  5  .±S  n0t  °Ut  °f  ran§e  wil1  Put  the  ’pen  down’  and 

f™tt0  he  callinS  Programme;  if  out  of  range  this  subroutine 

Coefficient  Pen"  ^  ^  Called  from  the  graphing  of  Extinction 

(b)  Subroutine  RERUN,  figure  37,  asks  for  a  run  number  and  then 

INSKOO  DAT66  •  Tfr0tn  diSC  int,°  mem°ry;  these  files  are 

INSKOO  DAT  no  smoke  ,  ISMK00.DAT  ’smoke’,  and  ICLKOO  DAT 

calcuiatW  (-°0.MT  is  a  ran  number  between  Si  and  5" 

(  background  is  not  necessary  here) . 

(c)  Subroutine  CLS,  figure  38,  clears  the  Visual  Display  Unit  (VDU) 
y  writing  twenty  five  lines.  Mode  switching  of  the  VDU  is  not 

satisfactory  since  several  VDU  types  may  be  used. 

(2)  MRL3 . FOR  listed  in  Appendix  III  is  the  filename  which  holds  the 

(a)  Subroutine  RAW(IFIL) ,  figure  39,  which  is  the  main  graph  drawing 

values 6 of  ’  f  thA  ®raph  -draWS  RdW  °ata’  a  sraph  of  absolute 

°  v  n.°  ,S1mokTf .  1  and  smoke  against  wavelength  position  of  the 
Circular  Variable  Filter  (CVF) .  (Comparison  of  this  with  previous 
s  s  can  point  to  system  misalignment  or  contamination.)  If  IFIL=1 
e  graph  drawn  is  for  Filter  Calibration  which  is  the  percentage 
-  smoke  x  100  s 


ratio 


'no  smoke'  VS  wavelength  position  of  the  CVF. 


H  IRIR  "  2  the  Mass  Extinction  Coefficient  is  drawn;  identical  to 
the  data  logger  programme  MRLV20. 


Subroutine  CLIN  is  very  similar  to  LOCLIN  except  the  limits  have 

cn an o-pH  UU  .r _  ,  .  .  .  .  . 


Cb) 

been  changed.  This  is  called  from  graphing  both 
filter  calibration'. 


raw  data'  and 


6.2  Data  logging  routines 


pie  main  data  logging  programme  MRLV20  is  formed  by  'linking'  four  editable 
files;  MRLV20.FOR,  MRLV21.F0R,  MRLV22 . FOR  and  READ43.MAC  flistings  are  Jn 
Appendices  I.IV.V  and  VI).  (There  is  a  variation  of  MRLV20  called  MRLF20 
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which  allows  greater  concentration  of  obscurant  from  -  100  g  to  -  30  kg  so 
fog  oils  may  be  tried  for  their  extinction  coefficients  (MRLV20's 
concentration  allowed  is  from  1  mg  to  100  g) . ) 

(1)  MRLV20.FOR,  figure  40,  comprises  the  MAIN  routine  which  consists  of 
a  'common  block'  of  RAM  (Random  Access  Memory),  fixed  data  block  for 
generating  run  numbers  for  dynamic  file  name  construction,  wavelength 
number  generation  equations,  initialising  routines.  It  reads  the 
current  run  number  stored  on  disc,  and  a  menu  which  calls  nine 
subroutines,  the  'no  smoke'  subroutine  twice. 

(a)  Subroutine  GRP,  figure  41,  is  the  actual  graphing  subroutine 
which  consists  of  the  'common  block',  it  assigns  the  plotter  handler 
to  'logical  unit  No  3',  and  then  prepares  for  plotting  -  such  as 
'paper,  pen,  plotter  ready?'.  (The  plotter  used  is  a  Houston 
Instrument  DMP-4  A4  size  single  pen  X-Y  plotter  having  a  step  size  of 
0.1  mm.)  After  defining  the  Origin  point  the  borders  are  drawn 
using  'tick'  markers  for  later  annotation.  Labelling  and  annotation 
is  followed  by  printing  the  appropriate  data,  and  finally  the  Mass 
Extinction  Coefficient  is  drawn.  The  initial  position  (X  and 
Y  coordinates)  is  determined  and  LOCLIN  is  called  to  position  and 
(normally)  to  put  the  'pen  down';  from  here  each  segment  is  drawn 
via  LOCLIN;  100/points  are  determined  per  segment  of  the  CVF,  and  if 
out  of  range  the  pen  is  'lifted'.  At  the  end  of  each  segment  the 
pen  is  lifted.  At  the  end  of  drawing  the  third  segment  the  raised 
pen  returns  to  the  origin.  Segments  numbers  2  and  3  were  installed 
in  reverse,  and  so  the  plotter  plots  these  two  from  right  to  left. 

(b)  Subroutine  LOCLIN  is  identical  to  the  one  described  in 
paragraph  6.1(a)  above. 

(2)  MRLV21.F0R,  figures  42  and  43  (listing  in  Appendix  V)  holds  the 
subroutines  RERUN  and  ERR0P. 

(a)  Subroutine  RERUN,  figure  42,  reloads  data  files  'nosmoke', 
'smoke',  and  calculations  stored  on  disc  back  into  memory;  and  is 
described  in  paragraph  6.1(b)  above. 

(b)  Subroutine  ERR0P(x,x,x) ,  figure  43,  an  error  detecting 
subroutine,  has  passed  to  it  two  variables  which  are  compared  to 
statistical  limits  which,  if  not  returned  to  zero,  give  a  confidence 
inference  of  90%  for  that  particular  plotted  reading. 

(3)  MRLV22.F0R,  figures  44  to  48  (listing  in  Appendix  VIII)  holds  the 
subroutine  CLS  which  clears  the  VDU  screen,  WEN  resets  the  run  number, 
NUR  alters  the  run  number,  NSK(IBAK)  does  'background'  and 

'no  smoke'  data  gathering,  SK  does  smoke  data  gathering,  CLC  calculates 
each  point,  TABUL  tabulates  results,  and  FIN  updates  run  number  before 
ending  the  programme  session. 

(a)  Subroutine  CLS  clears  the  VDU  screen  by  writing  25  blank  lines 
to  the  screen.  This  subroutine  is  identical  to  the  one  described  in 
paragraph  6.1(c)  above. 
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(b)  Subroutine  WEN*,  figure  44,  is  used  when  installing  a  new  data 
disc.  It  resets  the  IRUN  counter  to  '  1 1  .  This  subroutine  reduces 
the  problem  of  over-running  the  IRUN  numbers,  ie  it  keeps  the  run 
numbers  below  50. 

(c)  Subroutine  NUR*,  figure  45,  allows  a  single  key-stroke  to 
increment  the  run  number,  or  to  change  the  current  run  number  for  the 
purposes  of  reloading  a  set  of  data  files  from  disc,  or 
updating/upgrading  data,  eg  using  existing  data  to  patch  up  an 
incomplete  set  of  data  files.  (A  calculation  file  may  need  to  be 
copied  so  a  particular  run  number  has  the  three  files  required  for 
reloading) . 

(d)  Subroutine  NSK(IBAK),  figure  46,  is  used  to  obtain  a  background 
file,  and  a  no  smoke  file.  If  IBAK=0  a  ’no  smoke'  run  is  implied; 
but  if  IBAK— 1  a  background  run  is  implied.  The  main  difference 
between  these  two  'set-ups’  is  that  the  background  needs  the 
Nernst  Glower  blanked  off  so  any  hot  spots/reflections  etc  in  the 
Spectro-radiometer  and  Smoke  Chamber  itself  are  noted  and  then  can  be 
subtracted  from  the  'no  smoke'  and  'smoke'  runs  before  they  are 
manipulated  or  stored  on  disc. 

Looking  at  figure  46  shows  minimal  differences  except  for  'averaging' 
and  'storing';  the  main  difference  is  blanking  off  of  the 
Nernst  Glower  source,  to  get  a  background- reading  set  of  data,  the 
'numbers'  of  which  vary  with  the  CVF  wheel  position. 

Since  time-of-day  will  pinpoint  the  run  data  this  needs  to  be  put  in 
by  hand  because  the  'line  time  clock'  built  into  the  computer  has 
been  disabled,  so  timing  in  the  Macro  data  gathering  loop  can  be 
guaranteed. 

To  gain  more  reliable  data,  averaging  at  each  wavelength  is  required, 
but  since  memory  space  is  limited  a  maximum  of  7  revolutions,  of  the 
CVF,  can  be  stored  in  RAM  before  averaging.  Normally  3  to 
5  revolutions  would  be  used. 

When  the  parameters  are  set,  the  Macro  READ43 . MAC ,  figures  51  to  58, 
COLDAT  (collect  data)  is  called(ref . 1) .  On  returning  from  COLDAT, 
each  wavelength  position  is  averaged  before  being  stored  in  the 
respective  disc  file.  This  subroutine  NSK(IBAK)  returns  control  to 
the  menu. 

(e)  Subroutine  SK,  figure  47,  is  similar  to  NSK(IBAK)  in  that  it  is 
used  to  collect  and  store  smoke/obscurant  data.  It  asks  for  time- 
of-day,  smoke  type  description,  temperature  and  humidity  in  the 
Smoke  Chamber  itself,  and  then  the  number  of  revolutions  of  the  CVF 
required  (3  to  5  revolutions)  before  calling  the  data  collecting 
Macro  COLDAT.  Each  wavelength  is  averaged  and  then  background  is 
subtracted  before  being  stored  in  the  smoke  array  ISMK(I),  and 
finally  stored  onto  disc. 

(f)  Subroutine  CLC,  figure  48,  calculates  the  Mass  Extinction 
Coefficient  when  given  the  amount  of  material  that  was  used  during 
the  smoke  run  .  The  Mass  Extinction  Coefficient  parameters  are 


*  NEW  and  RUN  may  be  interpreted  by  the  computer  as  command  words 
which  could  make  it  crash  ,  so  the  simplest  method  is  to  reverse  the 
words  to  WEN  and  NUR  respectively. 
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tested  in  the  ERROP  subroutine  to  determine  if  the  reading  is  within 
the  90%  confidence  area,  if  not  it  will  be  forced  to  zero  where  the 
graphing  routine  will  detect  it  as  an  error  and  not  graph  that  point, 
but  leave  a  hole  in  the  graph  at  that  position.  The  ACALC(I)  data, 
date,  time,  concentration,  temperature  and  humidity  are  then  stored 
onto  disc  before  returning  to  the  menu. 

(g)  Subroutine  TABUL,  figure  49,  writes  all  the  data  to  the  VDU. 
It  includes  number  of  revolutions  used,  date,  time  of  'no  smoke'  and 
'smoke',  temperature,  humidity,  smoke  description,  concentration,  and 
the  three  files  of  data  ACALC(I),  INOSMK(I),  ISMK(I)  against 
ALAMDA(I). 

(This  part  of  the  programme  has  hard-copy  printout  'commented  out' 
since  there  is  no  readily  usable  printer  in  the  Smoke  Chamber 
building. ) 

The  programme  then  returns  to  the  Menu. 

(h)  Subroutine  FIN,  figure  50,  increments  the  IRUN  counter  which  is 
then  stored  onto  disc,  after  which  the  programme  session  (data 
logging)  ends . 

To  completely  describe  the  MRL  Smoke  Chamber  Instrumentation  and  Recording 
System  the  Macro  data  gathering  and  storing  subroutines  are  flow  diagrammed  in 
figures  51  to  58. 

A  breakdown  of  each  flow  diagram  is : - 
Figure  51  - 

The  main  controlling  Macro  CQLDAT,  with  parameters  passed,  and  called  from  the 
Fortran; 

Figure  52  - 

Subroutine  READ.  It  reads  and  stores  one  CVF  wheel  of  data; 

Figure  53  - 

PRMT  is  a  short  subroutine  which  displays  a  prompt  on  the  VDU; 

Figure  54  - 

Subroutine  ECD  is  the  'edge  change  detector'  which  detects  a  HEDS1000  Read 
Pulse ; 

Figure  55  - 

Subroutine  CNVT  performs  an  analogue  to  digital  conversion  when  required; 
Figure  56  - 

Subroutine  TRS  controls  the  data  gathering  for  100  pulses  or  readings  each 
segment; 

Figure  57  - 

Subroutine  TIM2  is  a  software  timer  which  is  used  in  taking  intersegment 
background  readings ; 

Figure  58  - 

Subroutine  TIM1  is  a  software  timer  for  1  ms  which  needs  to  be  adjusted  to 
suit  whatever  LSI  system  is  used.  Each  LSI  system  has  a  specific  data  clock 
rate  which  determines  the  instruction  timing. 
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DRCS  •  5  14  7530  TN  - 998 - 6616 
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APPENDIX  I 


READ43 . MAC 

.TITLE  READ43 . MAC 

9 


$ 

OLD 

.TITLE  RE AD4  >  MAC 

9 

SLINK 

y 

OB JS  OF 

MRLV20 .  FOR  .  MRLV21 .  FOR  . MRL V22 .FOR > READ43 . MAC 

. GLOBL 

COLDAT 

. MCALL 

.PRINT  .  .EXIT >  ♦  TTYOUT  . .TTYIN 

9 


start: 


MOV 

*  Q  .  R  2 

.INITIALLY  JAMS  R2  LOU 

MOV 

R  2  .  R 1 

JSR 

PC  » ECD 

5  FINDS  RISING  EDGE 

JSR 

PC. ECD 

. FINDS  FALLING  EDGE 

MOV 

0*CTR»R4 

CMP 

'  *300. .  @  *  C  T  R 

.LONG  CNT  IS  SYNC  PULSE 

BPL. 

START 

JSR 

PC . PRMT 

5 

BR 

START 

5 

CRL F3 : 

JSR 

PC. CRLF 

CRLF2 : 

JSR 

PC. CRLF 

crlf: 

TSTB 

0*177564 

BPL 

,  -4 

MOV 

*15 >0*177566 

TSTB 

0*177564 

BPL 

.  -4 

M  0  V 

*1-2  >0*177566 

RTS 

PC 

prmt: 

» PRINTS  PROMPT 

JSR 

PC. CRLF 

TSTB 

0*177564 

BPL 

.  -4 

MOV 

*45.0*1 77566 

.  % 

TSTB 

0*177564 

BPL 

.  -4 

MOV 

*40.0*177566 

>  SPACE 

A 

f 

RTS 

PC 

c  i  n  : 

.GETS  CHARACTER  IN  K/BOARD 

TSTB 

0*177560 

BPL 

.-4 

MO  VB 

0*177562. RO 

BIC 

*177600 . RO 

.WIPE  OFF  PARITY  BIT 

RTS 

PC 

co; 

.PUTS  CHAR  OUT  TO  VDU 

TSTB 

0*177564 

BPL 

.  -4 

MO  VB 

RO. 0*177566 

RTS 

PC 

5 


TIM1:  >  1  MILLISECOND  TIMER 
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MOV 

*60 .RO 

lpi  : 

COM 

R2 

COM 

R2 

DEC 

RO 

BNE 

LPI 

S 

RTS 

PC 

TIM25 

MOV 

*2220  . R5 

»  20  PULSES  D E L A Y  = 1 ♦ 17  SECS 

l.F'4  : 

JSR 

P  C  »  T I M 1 

DEC 

R5 

BNE 

LP4 

*» 

RTS 

PC 

cnvt: 

MOV 

*1  ? 0*177000 

5ADC  CONVERSION  SSTORE  DATA 

TSTB 

0*177000 

BPL 

.-4 

MOV 

0*177002.  -<R4)  + 

t 

MOV 

O  +  COUNT  t  <  R 4  )  + 

i  LLL 

INC 

0  *  C  0  U  N  T 

r  TEMP 

PATCH  TO 

OVERRIDE  ADC  NON 

EXISTENT  CARD 

5 

RTS 

PC 

$ 

TRS  : 

MOV 

*100  ,  • R5 

•TAKE  READINGS 
f TAKE  1  0  0  READINGS 

LP5J 

JSR 

PC. ECD 

•DETECTS  RISING  EDGE 

JSR 

PC . CNVT 

.ADC  AND  STORE  DATA 

JSR 

PC. ECD 

.DETECTS  FALLING  EDGE 

DEC 

R5 

BNE 

LP5 

*? 

RTS 

PC 

ecd: 

MOV 

*0.0*CTR 

.ZEROS  LOOP  COUNTER 

lp7: 

JSR 

PC.TIM1 

CMP 

*17000 . O+CTR 

.MOD  COUNT 

BMI 

LP9 

l.pi  i : 

MOV 

0*167774 . R 1 

’FILTER  WHEEL  INPUT 

BIC 

*77777. R1 

.BIT  15  USED  AS  INPUT 

CMP 

R 1  .  R2 

BEN 

LPIO 

MOV 

R1  .  R2 

RTS 

PC  • 

LP9 : 

DEC 

BR 

@*CTR 

LPI  1 

.MAINTAINS  CNTR  MAX  CNT 

lpio  : 

INC 

BR 

O  +  CTR 

L  P  7 
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5  BEFORE  THIS  PROGRAM  IS  USED  DFTSTR  INTO  R4 
5R4  THEN  CURRENT  STORE  LOCATION 


i READS 

read: 


L.P14 : 


LP.13 

9 


ONE  RESOLUTION  OF  FILTER  WHEEL  POINTS 


300  +  3  BACKGROUND 


MOV  *0?R2 

MOV  R 2  f  R  .1 

JSR  PCjECD 

JSR  F'CfECD 

CMP  *300.»e#CTR 

BPL.  READ 

MOV  ♦  3  f  R  3 

MOV  ♦  0  .*  R 2 

JSR  F'CfTRS 

JSR  PC»TIM2 

MOV  # 0 r ( R 4  )  + 

JSR  PCfCNVT 

MOV  *0»(R4>+ 

DEC  R  3 

BNE  L  P 1 4 

MOV  *0»(R4)+ 

MOV  ♦Of ( R 4  )  + 

RTS  PC 


» INITIALLY  JAMS  R2  LOW 

•FINDS  RISING  EDGE 
*  FINDS  FALLING  EDGE 
LONG  CNT  IS  SYNC  PULSE 

5  THREE  SEGMENTS  COUNTER  R3 

5  JAMS  R2  LOW  (SHOULD  ALREADY  BF  LOW) 
STAKE  100  DETECTOR  READINGS 
520  PULSES  DELAY  1,17  SECS 
•SEPARATOR 

> STORES  BACKGROUND  READING 
5  SEPARATOR 


5  SEPARATORS  DEPICTING  END  OF  REV 


coldat: 


•  MOV 

TST 
MOV 
MOV 
MOV 
MOV 
M  0  V 
MOV 

CL5  : 

MOV 

DEC 

BNE 

INC 

REV1  : 

DEC 
BEQ 
J  SR 
JSR 
BR 

RE VEND : 

RTS 


#0  .*  @*167770 
(  R  5  >  + 
(R5)+fRl 
!?(R5)  +  .R0 
■R 1  f  R4 
ROf ©♦REVS 
♦2200 , • RO 
♦OfGtCOUNT 

♦  Of  ( R  1  )  + 

RO 
CL  5 

©♦REVS 

©♦REVS 
RE VEND 
PC f READ 
PC  f  PRMT 
REV  1 

PC 


•READS  %  STORES  'ONE  REVOLUTION'  OF  DATA 
•TIMES  THE  NUMBER  OF  REVS  REQD. 

5  DISABLES  INTERRUPTS 

i DUMMY  (USED  TO  PASS  NUMBER  OF  PARAMETERS) 

•  T  E  M  F’ 

5  REVS 

•  F'  U  i  S  STORE  START  A  D  D  R  I N  R  4 

5  STORES  CURRENT  NUMBER  OF  REVS 

5  TEMPORARY  COUNTER  FOR  ADC  DATA  INPUTTING 

•  CLEARS  T F H p  STORE  A R E A 


•MARKS  END  OF  EACH  REVOLUTION 


count; 
ctr  : 

REVS  : 


,  BL.KW 


.BLKU 


1 

1 
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.  BL  KW  1 
.EVEN 


f 


END 


START 


n o n  n  o o n n  non  no 
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APPENDIX  II 
MRLGRP . FOR 

c  filename:-  mrlgrp. for 

C  OLD  FILE  NAME  FILENAME!-  MRLMAN.FOR 

C 

C  LINK  OB JS  OF  MRLGRP . FOR » MRL3 . FOR 

0 

0  MODIFIED  FOR  GRAPHING  ONLY 

C 

C  MRL  SMOKE  CHAMBER  PROGRAMME 

C  ACCEPTS  GOLAY  CELL  INPUT 

C  STORES  &  ANALIZES  DATA 

C 
C 

C  X 

COMMON  INQSMN (312) r ISMKI3.12) , ACALC(312) >  A  L  A  M  D  A  ( 103)  , 

#ALAMD<  103)  , AL AM  <  .105)  »RING(  10)  t  DTI  ,  DT2  ,  DT3  »  ANSKT  t  INSKR,  ASMKT,  ISMKR 
*t  TH,RH,CONC,  IBAKGD,  I  RUN  ,NRUN(50)  ?  I X  X  .•  I  Y  Y 
C 

LOGICAL *4  ANSKT  ,  A  S  M  K  T  »  R I N  G 
LOG  I  CAL *4  D  T 1  ,  D  T  2  »  D  T  3 

n 

o 

C.  DYNAMIC  GENERATION  OF  WAVELENGTH  NUMBERS 

C 

C  AMCEN 1=3.4155 

AMCEN 1 =3 . 5 1 25 
C  AMC£N2=6 .2020 

AMCEN2=6  *2762 
AMCEN3=11 .4972 
AMCEN3=11 ,3072 
AMCEN3=11 ,4381 
A  M  D  E  L 1 =  0 , 0200707 
AMDEL2=0 .0372490 
A  M  D  E  L  2  =  0, 0371. 1. 1  .1. 

AMDEL3=0 , 0659303 
AMDEL3-0 . 0566809 

A  L  A  M  D  A ( 1 )-(AMCENl~AMDELl / 2 . )-49< *AMDEL1 
ALAND  (  1  )  -•  <  AMCEN2  FAMDEL2/2  ,  >  +49  .*AMDEL2 
ALAM( 1 )-( AMCEN3+AMDEL3/2. ) +49 , *AMDEL3 
DO  530  1-29 100 

A  I  - 1 

AL AMDA ( I )=ALAMDA< 1 >+AMDELl*AI 
ALAMDt  I )=ALAMD< 1 )-AMDEL2*AI 
ALAMI I )=ALAM< 1 ) - A  M  D  E  L  3  *  A I 
530  CONTINUE 

END  OF  DYNAMIC  WAVELENGTH  GENERATION 


DATA  NRUN/'Ol ' >  ' 02  '  ,  '  03  '  t  '  04  '  ,  ' 05 ' ,  '  06  '  t  '  07  '  r  '  08  '  ,  ' 09 ' t  '  10  '  t 
+  '  1. 1  '  »  '  1  2  '  »  '  1  3  '  .  '  1 4  7  f  '  1 5  '  ,  '  1 6  '  »  '  1 7  '  ,  '  1 S  '  ,  '  1 9  '  t  '  20  '  , 

+  '21'j'22'»'23,,'24'>'25','26'»'27,?'28','29'»'30'i> 

+  '31' , '32' , '33' , '34' • '35' , '36' , '37' , '38' ,  '39' , '40' , 

+  '41' , '42' »  '43'  »  '44' » '45' »  '46' , '47' , '48' » '49' ,  '50'/ 


INITIALLIZE  COMMON  VARIABLES  TO  ZERO 


DO  3  1-1,10 


no  non 
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RING< I >=0 
3  CONTINUE 

C 

I NSKR=3 
ISMKR=3 
TH=00 , 0 
R  H  =  0  0  *  0 
CONC^OO , 0 
I B  A  K  G  D  =  0 


CALL  ASSIGN  (  2  ,  '  OY 1  J  J RUN  .  DAT  '  ,  0  ,  '  OLD  '  ) 
READ  (2,9)  I  RUN 
9  FORMAT  (12) 

CALL. CLOSE  (2) 


CALL  CLS 
TYPE  10 

10  FORMAT  <//.//) 

TYPE  13 

.13  FORMAT  (10X,'G  R  A  P  H  I  N  G  R  0  U  T  I  N  E  S'/) 

I RUN= 1 
80  TYPE  30 

30  FORMAT (  /  2  0  X  , ' M  E  N  U ' /20X •'=  =  =  =  =  =  -  =  =  -'/ ) 

IF  (  I R U N  »  G T  *  4 ? )  I R U N  =  1 
TYPE  35,. I  RUN 

35  FORMAT  (10X,  '  RUN  NUMBER  IS  #',12/) 

TYPE  40 

40  FORMAT  (5Xr'GIVE  A  DIGIT  BETWEEN  1  l  5'/) 

TYPE  30 

X 

50  FORMAT  (10X,1  -  LOAD  DISC-STORED  DATA  INTO  R  A  M'/lOX,  '2  =  GP 

*APH  EXTINCTION  COEFFICIENT') 

TYPE  60 

60  FORMAT  ( 10X,  '3  -  GRAPH  RAW  P  AT  A  '  /  .1  OX  ,  '  4  =  GRAPH  FILTER  CALIBRAT 

*  I  ON ' / 1  OX , ' 5  =  END  PROGRAMME '  ) 

READ  (5,*)  JA 

IF  ( JA.LT . 1 .OR. JA.GT  .5  )  GO  TO  80 
GO  TO  (100,105,110*  115,999)  J  A 
GOTO  SUBROUTINES 
105  I F I L  =  2 

CALL  RAW(IFIL) 

GO  TO  160 
100  CALL  RERUN 
GO  TO  160 
110  IFIL*0 

CALL  RAW  (I  FID 
GO  TO  160 

115  I F I L  =  1  - 

CALL  RAW ( IFIL  ) 

GO  TO  160 
160  CALL  CLS 

DO  190  I  =  1,10000 
190  CONTINUE 
GO  TO  SO 
999  CONTINUE 
TYPE  998 

998  FORMAT  (IX, 'END  OF  PROGRAMME') 
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C 

C 

C 

C 

C 

e 


C 


330 

335 

0 


300 
30  5 
310 


C 

C 

r' 

U 

C 

c 


c 


c 

c 


c 


705 

710 


711 

C 


715 

716 


C 


720 


END 


LOCATES  START  OF  LINE  AND  PUTS  PEN  DOWN 
SUBROUTINE  LOCLIN 

COMMON  INQSMK ( 312 ) . ISMK (312) . ACALC(312) *  A  L A  M  D  A ( 103)  . 

*ALAMD< 103) . ALAM( 105) >RING( 10) .  DTI . DT2.  DT3. ANSKT  . INSKR .  AS MKT . I SMI'R 
*.TH.RH.CONC. I B AKGD » I  RUN . NR UN ( 50 ) . I XX . I Y Y 

IF  (IYY.GT *1200*0 R*IYY.LE*0)  GO  TO  335 
WRITE  (  3  »  3  0  0  )  I X  X  . 1 Y  Y 
WRITE  (3.305) 

PEN  DOWN. 

RETURN 

WRITE  (3.310) 

PEN  UP 
GO  TO  330 

F  0  R  M  A  T  (  A .16.  .'.16.'  ) 

FORMAT  (  '  D  '■  ) 

FORMAT ( ' U  ') 

END 


RE-RUN  OF  EXISTING  DATA  STORED  ON  DISC 
SUBROUTINE  RERUN 

X 

COMMON  I  NOSH i<  (  3 1 2  >  .  I SMK  (  3  I.  2  )  ;  ACALC  (312)  .  ALA M D  A  (103). 

* ALAMD ( 103) . ALAM ( 105) . RING ( 10) .DTI . DT2 . DT3 . ANSKT  >  INSKR  .  ASMKT . I SMKR 
*  .  TH . RH . CONC ? IBAKGD. I  RUN . NR UN ( 50) . I XX . I YY 

I  N  T  E  G  E  R  *  2  D  I N  K  (  7  )  .  D  I  S  K  (  7  )  .  D I C  K  (  7  ) 


D  A  T  A  D I N  K  /  D  T  .  1  *  ,  .  '  I  N  '  .  '  S  K  '  .  0  0  '  .  '  <  D  *  ‘  fa  j  '  / 

DATA  D 1 3  K  /  '  D  Y  '  .  '  1  5  '  *  '  1 5  '  .  '  H  K  '  .  '  0  0  '  .  '  *  D  '  A  T  '  / 

D  A  T  A  DIC  K /  D  Y ' >  ' 1 »  ' »  ' I C  '  .  '  L  N  '  .  '  0 0  '  .  '  *  D  '  *  AT'  / 

CALL  CLS 
TYPE  710 

FORMAT  (10X.  'RE-RUN  OF  PREVIOUS  DISC-STORED  DATA '//5X»  '  GIVE 
*R  U  N  NUMBER  YOU  WANT  TO  SEE') 

READ  (5.711)  IRUN 
FORMAT  (12) 


DINK(5)=NRUN( IRUN) 

CALL  ASSIGN  ( 2 . D I NK . 1 4 . ' OLD '  ) 

READ  (2.715)  ( I NOSHK ( I ) . I = 1 . 3 1 2 ) 
FORMAT  ( 63 ( 5  < 1 1 0 . 2X ) ? / )  ) 

READ  (2.716)  DTI , DT2 » DT3 » ANSKT » INSKR 
FORMAT  (A2.A3.A2.5X.A4.5X. 12) 

CALL  CLOSE  (2) 


D I S  K ( 5 ) =  N  R  U  N ( IRUN) 

CALL  ASSIGN  ( 2 . D I SK . 1 4 , ' OLD '  ) 
READ  (2.720)  (ISM  K ( I ) . 1-1.312) 
FORMAT  (63(5(1 10 . 2X) , /)  ) 


n  n  n 
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READ  (2,721)  DTI , DT2 » DT3 » ASMKT , I SMKR , ( RING ( I ) , 1=1 , 10 ) 

721  FORMAT  ( A2 , A3 , A2 , 5X , A4 , 5X , 12 , 1 0A4 ) 

CALL  CLOSE  (2) 

C 

DICK(5)=NRUN( IRUN) 

CALL  ASSIGN  < 2 , D ICK , 1 4 , ' OLD  '  ) 

READ  (2,725)  ( ACALC ( I ) , I  =  1 , 3 12 ) 

725  FORMAT  (63(5(F6.3,2X) ,/)  ) 

READ  (2,726)  DT 1 , DT2 , DT3 , ASMKT , CONC , TH , RH , I BAKGD 

726  FORMAT  (A2»A3,A2,5X,A4,5X,F6*4,5X,F6«2,5X,F6*2,5X,I6) 

CALL  CLOSE  (2) 

C 

TYPE  730, IRUN, DTI ,DT2,DT3, (RING (I), 1=1, 10) ,CONC 
730  FORMAT  (//5X,'RUN  *', 12 , 3X ,' DATE l - '  , A2 , A3 , A2/5X  ,'  SMOKE  TYPE,-'  , 

♦10A4/5X, 'CONCENTRATION  ',F6.4,'  GMS/CUB . METRE ' // 1  OX ,  ' I S  THIS  THE 
♦CORRECT  RUN  ?  Y/N') 

READ  (5,735)  COR 
735  FORMAT  (Al) 

IF  ( COR , NE  .  ' Y' )  GO  TO  705 

RETURN 

END 


SUBROUTINE  CL.S 

C  CLEARS  0  D  U  SCREEN 

TYPE  5 

5  FORMAT  (//////////.///////////////) 

RETURN 

END 
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APPENDIX  III 


C 

C 

C 

C 

c 

c 

c 


c 

c 

c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


c 

c 

c 

c 

c 


501 

510 


300 


305 

512 

C 

500 

504 

505 


520 


525 


540 

545 


C 

C 

C 


MRL3 . FOR 

FILENAME:-  MRL3.FQR 
LINK  OB JS  OF  MRLGRP.F0RfMRL3.F0R 

MRL  SMOKE  CHAMBER  PROGRAMME 
ACCEPTS  GOLAY  CELL  INPUT 
STORES  8  ANALIZES  DATA 


X 


DRAWS  A  GRAPH  OF  RAW  DATA  IFIL=0 
DRAWS  FILTER  CALIBRATION  IF  1 1  =3 
DRAWS  EXTINCTION  COEFFICIENT  IFJL -2 


SUBROUTINE  RAW(IFIL) 


COMMON  I N  0  S  M K ( 3 1 2 ) , ISMK(312) , ACAL( 
♦ALAMD ( 103 ) , ALAM ( 105 ) , RING ( 10 ) , DTI 
*’THfRHf CQNC , IBAKGD  , I  RUN , NRUN ( 50 ) 


f  DT 
,  IXX 


(312) >  A  L  A  M  0  A  ( 103) , 

’  DT3 , ANSKT  , INSKR »  ASMKT , ISMKR 
IYY 


CALL  A3SIGN(3,  'PL:  ' >0 ,  'NEW' ,  '  N  C  '  } 

CALL  A  S  3 1 G  N  (  3  , '  T  T :  '  , 0  »  'NEW' ,  ' N  C ' ) 

CALL  ASSIGN (3,  'LP:  '  ,0,  'NEW' >  'NC  '  > 

CALL  CL.S 

FORMAT  ( / / / / 2  0  X  ?  ' F  I  L  T  E  R  CALIBRATIO  N  ' /// 1  OX >  '  PAPE 
Y  R  READY  Y / N  ?  '  / ) 

IF  (IFIL.EQ.O)  GO  TO  512 
IF  (IFIL.EQ.2)  GO  TO  300 
TYPE  501 
GO  TO  504 
TYPE  305 
GO  TO  504 

FORMAT  <////10Xf'G  R  A  P  H  I  N  G  EXTICTION  C  0  E  F  F' 
*///10X f  '  PAPER  READY  Y/N  ?'/) 

TYPE  500 

x 

FORMAT  (//// 20X »  ' G  R  A  P  H  I  N  G  R  A  W  D  A  T  A  '  / / / 1 0  X  * ' PAPER 
♦READY  Y/N  ?'/) 

READ  (5,505)  PAP 
FORMAT  <  A 1 ) 

IF  ( PAP . NE. ' Y '  )  GO  TO  510 
TYPE  520 

FORMAT  (lOXf'PEN  READY  Y/N  ?'/) 

READ  (5,525)  PEN 
FORMAT  (Al) 

IF  (PEN.NE.'Y')  GO  TO  510 
TYPE  540 

FORMAT  <10X, 'PLOTTER  ACTIVE  Y/N  ?  ' / ) 

READ  (5,543)  PLT 
FORMAT  (Al) 

IF  (PLT . N E .  'Y')  GOTO  530 


DRAW,  PLOT,  ETC 
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C 

C 

C 

650 

C 

C 

C 


590 

595 


C 

C 

C 


605 

600 

C 


C 


615 

610 

C 


C 


625 

620 


C 


635 

630 


C 

c 


c 

660 


C 

c 


ACTIVATE  PLOTTER 
WRITE  (3*650) 

PLOTTER  ACTIVE*  HOME (  EOT . L/H  CNR ) *  POS I T ION  2  C M  UP  AND 
ALONG*ORIGIN  AT  THAT  POINT 
FORMAT  <  '  H  A100 *  200  0  '  > 


GO  TO  888 

IF  (IFIL.EQ.2)  60  TO  310 
FORMAT! 'R ' *  15* ' * ' *  15*  '  ') 

FORMAT ( '  M20  ') 

595  =  TICK  FOR  X-AXES  AND  Y-AXES 


DRAW  BORDER  OF  GRAPH 
WRITE  (3*565) 

DO  600  1=1*15 
WRITE  (3*595) 

WRITE  (3*565) 

WRITE  (3*605) 

FORMAT!  ' R 150*0  '  > 
CONTINUE 

DO  610  1=1*10 
WRITE  (3*595) 

WRITE  (3*565) 

WRITE  (3*615) 

FORMAT ( 'RO* 100  ') 
CONTINUE 

DO  620  1=1*15 
WRITE  (3*595) 

WRITE  (3*565) 

WRITE  (3*625; 

FORMAT! 'R- 150*0  ') 
CONTINUE 


RAW  DATA  X  FILTER  CAL 


DO  630  1=1*10 
WRITE  (3*595) 

WRITE  (3*565) 

WRITE  (3*635) 

FORMAT!  'RO* -100  ') 

CONTINUE 

WRITE  (3*575) 

END  OF  BORDER 
L  A  B  =  0 
IXX=-100 
I Y  Y  •--  -  2  5 

WRITE  (3*560)  I X  X  *  I Y  Y 
DO  660  1-1*10 
WRITE  (3*655)  LAB 
LAB  =  l.AB+l 
WRITE  (3*615X 
WRITE  (3*710) 

STEP  VERTICAL 
CONTINUE 

WRITE  (3*655)  LAB 
GO  TO  215 

DRAW  EXTINCTION  BORDER 


on  on  o  no 
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310  WRITE  (3*565) 

DO  315  1=1*15 
WRITE  (3*595) 

C  WRITE  (3*565) 

WRITE  (3*605) 

315  CONTINUE 
C 

DO  316  1=1*6 
WRITE  (3*595) 

C  WRITE  (3*565) 

WRITE  (3*317) 

317  FORMAT (  '  RO  *  200  '  ) 

316  CONTINUE 
C 

DO  3 IS  1  =  1*15 
WRITE  (3*595) 

C  WRITE  (3*565) 

WRITE  (3*625) 

318  CONTINUE 
C 

DO  319  1=1*6 
WRITE  (3*595) 

C  WRITE  (3*565) 

WRITE  (3*320) 

320  FORMAT ( ' RO  * -200  ') 

319  CONTINUE 
WRITE  (3*575) 

C 

C  END  OF  BORDER 

L  A  B  =  0 
I X  X  =  -  8  0 
I Y  Y  =  -  2  5 

WRITE  (3*560)  I XX* IVY 
DO  210  1=1*7 
WRITE  (3*655)  LAS 
LAB=LAB+1 
WRITE  (3*317) 

WRITE  (3*710) 

210  CONTINUE 

WRITE  LABELLING 

215  L AB=0 

I XX=-60 
I Y  Y  =  -  5  0 

WRITE  (3*560)  IXX*IYY 
POSITION  ABSOLUTE 
DO  645  1=1*15 
WRITE  (3*653)  LAB 
655  FORMAT ('S12  '  *  I  2  *  '  _  '  ) 

HEIGHT  =  2. 5mm* ROTA TION  =  NORMAL 

LAB=LAB+1 

WRITE  (3*705) 

MOVE  TO  NEXT  LOCATION 
645  CONTINUE 

WRITE  (3*655)  LAB 
705  FORMAT ( '  R 1 0  5  *  0  ') 

710  FORMAT ( ' R-50  *  0  '  ) 

ANNOTATION 

C 


noon  noon 
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665 

C 


220 

670 

C 

230 


C 

675 


680 


685 


420 


235 

240 

410 

690 


IXX=650 

IYY=-120 

WRITE  (3.560)  IXX.IYY 
WRITE  (3.665) 

FORMAT (  '  S12  WAVELENGTH  (uM)_') 

HEIGHT-3  mm, ROTATION -  NORMAL. 

IXX--100 
I Y  Y-200 

WRITE  (3.560)  IXX.IYY 
IF  (IFIL.NE.2)  GO  TO  220 
WRITE  (3.225) 

FORMAT ('S42  EXTINCTION  COEFFICIENT  (  S  Q  M  /  G  M )  _  '  ) 

GO  TO  230 
WRITE  (3,670) 

FORMAT  (  '  S  4  2  NORMALLIZED  DATA  (N/1.())_') 

HEIGHT-3 mm, ROTATION -270" 

IXX-0 

IYY=1230 

WRITE  (3,560)  IXX.IYY 
WRITE  (3,675)  TH.RH 

FORMAT ( '312  TEMPI-' ,F6,2. 'C  REL . HUMIDITY J - ' ,F6.2, ) 

IXX-450 

IYY=1350 

WRITE  (3,560)  IXX.IYY 

WRITE  (3,680)  DTI , DT2.DT3.A3MKT , IRUN 

FORMAT ( '312  DATE:-  ' , A2 , ' / ' , A3 , ' / ' A2 , '  TIME:~',A4,' 

##' , 12, '_' ) 

I  XX-200 
I Y  Y  =  1 2  9  0 

WRITE  (3,560)  IXX.IYY 

WRITE  (3,685)  (RING (I), 1-1,10) 

FORMAT  ('312  SMOKE  TYPE  ♦  -  1 0A4  ,'...' ) 

I XX- 150 
IYY-1450 

WRITE  (3,560)  IXX,IYY 

FORMAT  ('313  FILTER  CALIBRATION-' ) 

IF  ( IFIL  ,EQ.O)  GO  TO  410 
IF  (IFIL.EQ.2)  GO  TO  235 
WRITE  (3,420) 

GO  TO  470 
WRITE  (3,240) 

FORMAT ('S13  MRL  SMOKE  CHAMBER  MEASUREMENTS-') 

GO  TO  470 
WRITE  (3,690) 

FORMAT ('S13  MRL  SMOKE  CHAMBER  RAW  DATA-') 


X 


RUN 


END  OF  ANNOTATION 
470  FIL-FIL 

560  FORMAT ('A', 15, ',',15,'  ') 

PLOTTING  ACTUAL  VALUES  FROM  DATA 


720  FORMAT (  ' M2 1  ') 

565  FORMAT ( ' D  ') 

C  PEN  DOWN 

575  FORMAT (  '  U  ') 


CJ  u 
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C  PEN  UP 

WRITE  (3,575) 

C 

IF  (IFIL.EQ.O)  GO  TO  888 
IF  ( IF IL . EQ . 2 )  GO  TO  350 

*  *  *  *  *  F I L  T  E  R  CALIBRATION**#*# 

I X  X  =  I N  T ( (ALAND A(3)*150)+0.5) 

AN0S  =  FL0AT(IN0SMK(3)  ) 

IF  (ANOS.EQ.O)  I Y Y=0 
IF  (ANOS.EQ.O)  GO  TO  400 
IYY  =  INT(FLQAT(ISMK(3))/AN0S*1000  +  0<.  5) 

400  CALL  CLIN 

DO  440  1=3,100 
I X  X  =  I N  T ( ( ALAMDA( I )*150 ) +0 . 5) 

ANOS  =  FLQAT( INOSMK (  I  )  ) 

IF  (ANOS.EQ.O)  I Y  Y  =  0 

IF  (ANOS.EQ.O)  GO  TO  401 

I Y  Y  =  I N  T  (  F  L  0  A  T  (  1 5  M  K  ( I )  )  /  A  N  0  S  *  1  0  0  0  +  0 . 5  ) 

401  CALL  CLIN 

440  CONTINUE 

WRITE  (3,575) 

C 

I  XX=  I  NT  (  (  AL.  AMD  (  3  )  #  1 50  )  +0 . 5  ) 

ANOS=FLOAT( INOSMK( 106) ) 

IF  (ANOS.EQ.O)  I Y  Y  =  0 

IF  (ANOS.EQ.O)  GO  TO  402 

IYY=INT  (FLOAT  (  ISMK<  106)  ) /ANGS*1()00  +  0 . 5  ) 

402  CALL  CLIN 

DO  450  1=106,202 

IXX=INT( ( ALAND < 1-103) *150) +0.5) 

C  IXX=INT ( (ALAND ( 1-105) *150) +0.5) 

A  N  0  S  =  F  L  0  A  T  ( I N  0  S  M  K  ( I  )  ) 

IF  (ANOS.EQ.O)  I Y  Y  =  0 

IF  (ANOS.EQ.O)  GO  TO  403 

I Y Y=  I  NT (FLOAT ( I  SMK ( I )  ) / ANOS* 1 000  +  0 . 5 ) 

403  CALL  CLIN 

450  CONTINUE 

WRITE  (3,575) 

C 

C 

I XX  =  I  NT ( <ALAM(3)#150)+0,5) 

AN03=FL0AT ( IN0SMK(209) ) 

IF  (ANOS.EQ.O)  I Y  Y  =  0 

IF  (ANOS.EQ.O)  GO  TO  404 

IYY=INT(FL0AT(ISMK(209))/AN0S*1000+0,5) 

404  CALL  CLIN 

DO  460  1=209,305 

IXX=INT ( (  A  L  A  M  ( I~206)*150)+0,5) 

C  I X  X  =  I N  T ( (ALAM( I~20S)*150)+0.5) 

ANOS=FLOAT ( INOSMK(I) ) 

IF  (ANOS.EQ.O)  I Y  Y  =  0 

IF  (ANOS.EQ.O)  GO  TO  405 

IYY=INT (FLOAT ( ISMK( I ) ) /ANOS* 1 000+0 , 5 ) 

405  CALL  CLIN 

460  CONTINUE 

WRITE  (3,575) 

GO  TO  430 
C 

C  *  *  *  *  *  N  0  S  M  0  K  E  THEN  SHORE  PLOT****' 


nn  n  nnon  n  no 
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888  I XX- 1 NT  < <ALAMDA(3>*150)+0.5) 

I Y  Y  =  I N  T (FLOAT ( IN0SHK( 3) )  /  2  +  0 . 5  ) 
CALL  CLIN 
DO  570  1=3,100 
I XX= INT ( <ALAMDA(I)*150)+0*5> 
IYY=INT ( FLOAT ( INOSMK ( I ) )/2+0.5> 
CALL  CLIN 
570  CONTINUE 

WRITE  (3,575) 


IXX=INT ( <ALAMD(3)*150)+0*5> 

I YY=INT ( FLOAT ( INOSMK ( 106) )/2+0*5> 
CALL  CLIN 
DO  580  1=106 ,  202 
I XX  =  INK  < ALAMO < I -103) *150) +0.5) 
IYY  =  .INT  (FLOAT  <  IN0SMK(  I )  )/2+0*5) 
CALL  CLIN 
580  CONTINUE 

WRITE  <3,575) 

IXX=INT< (ALAM(3)*150)+0*5) 
IYY=INT< FLOAT < INOSMK (209) )/2+0.5) 
CALL  CLIN 
DO  585  1=209,305 
IXX=INT<  <AL.  AM  (1-206)  *150)  +0.5) 
IYY=INT < FLOAT < INOSMK < I ) )/2  +  0*5  ) 
CALL  CLIN 
535  CONTINUE 

WRITE  (3,575) 

PEN  UP  -  AT  END  OF  PLOTTING 


88  I X  X  = I N  T ( ( ALAMDA(3)*150)+0*5) 
IYY=INT (FLOAT ( ISMK(3) )/2+0*5) 
CALL  CLIN 
DO  70  1=3,100 

IXX  =  INT(  (ALAMDA(  I)  *150+0*5) 

I Y  Y  = I N  T (FLOAT ( I S  M  K ( I ) ) / 2  +  0 . 5  > 
CALL  CLIN 
70  CONTINUE 

WRITE  (3,575) 

IXX  =  INT ( ( A  L  A  M  D ( 3 ) *  1 5  0 ) +  0 . 5 ) 
IYY=INT (FLOAT ( I S  M  K  < 106) )/2  +  0.5) 
CALL  CLIN 
DO  80  1=106,202 

I XX  =  I  NT ( ( ALAMO ( I - 1 03 ) *  1 50 ) +0 . 5 ) 
IYY=INT ( FLOAT ( I S  M  K ( I ) ) / 2  +  0 ♦ 5) 
CALL  CLIN 
SO  CONTINUE 

WRITE  (3,575) 


I XX= I  NT ( (ALAM(3)*150)+0.5) 

IYY=INT (FLOAT ( ISMK(20?) )/2+0.5) 

CALL  CLIN 

DO  35  1=209,305 

IXX  =  INT( <ALAM< I -206 ) * 1 50 ) +0  *  5 ) 


nnnn  nnnnoonn  n  n_  no 
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IYY=INT (FLOAT (  ISMK( I )  ) / 2  +  0 . 5 ) 

CALL  CLIN 
S5  CONTINUE 

GO  TO  430 

*****EXTINCTI0N  PLOT******* 
50  I XX= I  NT ( (ALAMDA(3)*150)+0»5) 

I YY= I  NT ( (ACALC(3)*200>+0.5) 

CALL  LOCLIN 

DO  355  I =3  * 100 

IXX-INT  < (ALAMO A( I )  *  1 5  0  )  +  0  ♦  5  ) 

I Y  Y  =  I N  T  ( (ACALC(I)*200)+0*5) 
CALL  LOCLIN 
55  CONTINUE 

WRITE  (3*575) 


I  XX=  I  NT  (  (ALAMD(3)*1. 50)+0»5) 
IYY-INT ( (ACALC(106)*200)+0*5) 
CALL  LOCLIN 
DO  360  1=106*202 
IXX=INT( (ALAMO ( 1-103) *150) +0.5) 
IYY=INT ( ( AC ALC ( I ) *200 ) +0  *  5 ) 

CALL  LOCLIN 
CONTINUE 
WRITE  (3*575) 


I X  X  =  I N  T  (  (ALAM(3)*.150)+0.5) 

I Y  Y  = I N  T ( (ACALC(20? ) *200) +0.5) 

CALL  LOCLIN 

DO  365  1=209*305 

IXX=INT( (ALAM( 1-206) *150) +0,5) 

I Y Y= I  NT ( ( ACALC( I ) *200) +0.5) 

CALL  LOCLIN 

CONTINUE 


WRITE  (3*575) 

PEN  UP  -  AT  END  OF  PLOTTING 


WRITE  (3*640) 

FORMAT  (  "H  (?  '  ) 

INACTIVATE  PLOTTER 
700  FORMAT( ' A' *  16* ' * ' *  16* '  ') 

CALL  CLOSE  (3) 


TYPE  550 

550  FORMAT  (///10X*'END  OF  GRAPHING '* 10X *'  PRESS  RETURN  KEY'/) 
READ  (5*555)  DUD 
555  FORMAT  (Al) 

RETURN 

END 


LOCATES  START  OF  LINE  AND  PUTS  PEN  DOWN 
SUBROUTINE  CLIN 


COMMON  IN0SMKC312) *ISMK(312) *ACALC(312) *ALAMDA( 103)  * 


X 


o  n 
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*ALAMD<  103)  »  ALAM(  105)  »RING(  10)  »DT1 »  DT2  >  DT3  »  ANSKT  ,  INSKR*  ASMKT  » ISMKR 
*, THf RH>C0NC> IBAKGD* I  RUN, NRUN(50) *IXX* IYY 
C 

IF  (IYY.GT.IOOO.OR.IYY.LE.O)  GO  TO  35 
WRITE  ( 3 »  30 )  I XX  *  I Y Y 
WRITE  (3.5) 

C  PEN  DOWN 

33  RETURN 

35  WRITE  (3*10) 

C  PEN  UP 

GO  TO  33 

30  F ORMAT ( ' A ' >  1 6  *  '  >  '  *  1 6 ,  '  '  ) 

5  FORMAT ( ' D  '  ) 

10  FORMAT (  '  U  ') 

END 
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APPENDIX  IV 


MRLV20 . FOR 


C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 


c 


c 

c 


c 

c 

c 


FILENAME:-  MRLV20.F0R 

OLD  filename:-  mrlman.for 

LINK  OB JS  OF  MRLy2O.FOR.MRLV21.F0R.MRLV22.FOR.READ43.MAC 

MRL  SMOKE  CHAMBER  PROGRAMME 
ACCEPTS  GOLAY  CELL  INPUT 
STORES  8  ANALIZES  DATA 

MODIFIED  TO  USE  312  BACKGROUND  READINGS 


X 

COMMON  I  TEMPI  2200) . TN0SMK(312) . ISMKI312 > . ACALC <  3 1 2 ) , A  L  A  M  D  A  ( 103) ? 
* AL AMD (  103)  *  AL AM (  105)  >  R I NG  <  10)  » DT 1  . DT2 . DT3  »  ANSKT  ■>  INSKR . ASMKT » ISMKR 
*  .  TH  .  RH  .  CONC  .  IBAKGD  .  I  RIJN  .  NR  UN  (  50  )  *  TXX  .  I YY  . 1  BACK  (312) 


LOGICAL *4  ANSKT > ASMKT  .RING 
LOGICAL *4  D  T 1 »  0  T  2  .  D T 3 


DATA  NRUN/'Ol ' * ' 02 ' . ' 03 ' .  '  04  '  >  '  05  ‘  ♦ / 06 
+  ' 11 ' . '12' . ' 13' . ' 14  ;  '15' .  '16 


07' . '08' . '09 
17' f '  18'  .  '  19 


10'  , 
20'  . 


+ 

+ 

+ 


'21 ' * '22' .  '23'  .  24  ' .  '25'  »  '26' .  '27' .  '28' . '29' . '30'  . 
' 31 ' . ' 32 • , ' 33 ' . ' 34 ' .  ' 35 ' . ' 36 ' . ' 37 ' . ' 38 ' . ' 39 ' . ' 40 ' , 
'41 ' .  42  '  .  '43' *  '44' . '45' . '46' . '47' » '48' > '49' .  'GO'/ 


DYNAMIC  GENERATION  OF  WAVELENGTH  NUMBERS 


C 

AMCEN1=3.5125 
A  M  C  E  N  2  =  6 . 2  7  6  2 
A  M  C  E N 3  =  1 1 .4381 
AnDEL 1=0 .0200707 
AMDEL.2  =  0, 0372490 
A  M  D  E  L  3  =  0 , 0659303 
C 

ALAMDA ( 1 )=( AMCEN1-AMDEL1/2. ) -49 . #AMDEL 1 
AL AMD  < 1 >  =  < AMCEN2  +  AMDEL2/2. )+49.*AMDEL2 
AL AM ( 1 )  =  ( AMCEN3  +  AMDEL3/2.  ) +49 . *AMDFL3 
DO  200  1-2*100 
A 1  =  I 

AL.AMDA(  I  )=ALAMDA<  1  )  +AMDEL 1  *A I 
ALAND ( I )=ALAMD(  1  >-AMDEL2*AI 
A L  A M  (  I  )  =  A L  A M  (  1  ) -AMDEL3*A I 
200  CONTINUE 
C 

C  END  OF  DYNAMIC  WAVELENGTH  GENERATION 

C  =  =  =  =  =  =  =  =  =:=:  =  =  =  =  =  =  =:  =  =  ^  =  =  =  =  =  =  =  =  =  =  =  ==  =  =  = 

p 

C  INITIALLIZE  COMMON  VARIABLES  TO  ZERO 

P 

w 

J  =  1 

00  1  1-1*2200 
J  =  J  +  1 

I  TEMP ( I )=J 
1  CONTINUE 

C 


DO  3  1=1.10 
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3 

C 


C 

C 

C 


9 

C 

C 

10 

.1.3 

18 

20 

80 

30 

35 

40 

50 

60 

70 


77 

C 

95 

100 

110 


ft  I NG ( I ) =0 
CONTINUE 

I N  S  K  ft  =  3 
I S  M  K  ft  =  3 
T  H  =  0  0 . 0 
ft  H  -  0  0  *  0 
C0NC=00 . 0 
I  B  A  K  G  D  =  0 


CALL  ASSIGN  ( 2 .  ' D Y 1 i I  RUN . D AT  '  , 0 . ' OLD ' > 
READ  (2;?)  I  RUN 
FORMAT  (12) 

CALL  CLOSE  (2) 


CALL  CL.3 
TYPE  10 
FORMAT  <////) 

TYPE  13 

FORMAT  (lOXi'S  MOKE  CHAMRER  MEASUREMENTS 
*'//> 

TYPE  18 

FORMAT  (IOXf'GIUE  DATE  EG  02N0V82'/) 

READ  (5.20)  D T 1 »  D T 2  f  D T 3 
FORMAT  (  A  .2  ?  A  3  .  A  2  ) 

TYPE  30 

FORMAT (  /  /  /  /  2  0  X  .  ' M  E  N  U ' /20X .  ' =  =  =====  =  =  = ' / ) 

IF  (IRUN.GT.49)  IRLJN=1 
TYPE  35 » I  RUN 

FORMAT  <10X.'  RUN  NUMBER  IS 
TYPE  40 

FORMAT  ( 5X  * ' G I UE  A  DIGIT  BETWEEN  1  &  10') 

TYPE  50 

FORMAT  (lOX.'l  =  TAKE  BACKGROUND  READINGS'/! OX. '2  =  NO  SMOKE  MEAS 
*UREMENT'/10X> '3  =  SMOKE  MEASUREMENT') 

TYPE  60 

FORMAT  (10Xf'4  =  CALCULATE  DATA ' / 10X , ' 5  =  GRAPH  RESULTS') 

TYPE  70 

FORMAT  (  1  OX  >  '6  -  TABULATE  RESULTS '/ 1  OX *' 7  -  NEW  DATA  DISKETTE  HAS 
*  BEEN  INSTALLED') 

TYPE  75 

FORMAT  (10X.'8  =  NEW  RUN  NUMBER  REQUIRED ' /10X »' 9  =  CLOSE  DOWN  CO 
*MPUTER  FOR  THE  DAY') 

TYPE  77 

FORMAT  (10X»'10=  RE-RUN  OF  PREVIOUS  DISC-STORED  DATA') 

READ  (5.*)  JA 

IF  ( JA . LT . 1 . OR . JA . GT  .  10)  GO  TO  80 
GO  TO  <95r 100. 110. 120. 130.140. 143. 145f 150. 155)  JA 
GOTO  SUBROUTINES 
1 .8  A  K  =  1 

CALL  NSK(IBAK) 

GO  TO  160 
I  B  A  K  =  0 

CALL  NSK(IBAK) 

GO  TO  160 
CALL  SK 
GO  TO  160 


in  m  in  in  in  to  in  in 

U  U  O  U  U  (J  G  U  O  U 
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120  CALL  CLC 
GO  TO  160 
1.30  CALL  GRP 
GO  TO  160 
140  CALL  TABUL 
GO  TO  160 
143  CALL  WEN 
GO  TO  160 
145  CALL  NUR 
GO  TO  160 
150  CALL  FIN 
GO  TO  999 
155  CALL  RERUN 
GO  TO  160 
160  CALL  CLS 

00  190  I  =  1,10000 
190  CONTINUE 
GO  TO  80 
999  CONTINUE 
TYPE  998 

998  FORMAT  (IX, 'END  OF  PROGRAMME') 

END 

C 

C 

C 

c 

c 

C  DRAWS  A  GRAPH  OF  EXTINCTION  COEFFICIENT  V  WAVELENGTH 

C 

SUBROUTINE  GRP 


COMMON  I  TEMP  (2200)  ,  I NOSMK  (  3 1 2  )  ,  I SMK  <  31 2  )  ,  AC AL..C  (  3 1 2  )  ,  ALAMDA  ( 103  )  , 
*ALAMD( 103) , ALAM( 105) , RING < 10 ) , DTI , DT2 , DT3 , ANSKT , INSKR  ,  ASMKT , ISMKR 
#,TH,RH,CONC, IBAKGD, IRUN,NRUN(50) , I XX , I YY , J B ACK ( 3 1 2 ) 


CALL  ASS  I GN ( 3 ,  ' PL t ' , 0 , ' NEW ' , ' NC ' ) 

CALL  A  S  S I G  N ( 3 , ' T  T  J ' ,0,  'NEW' , 'NC' ) 

CALL  ASSIGN (3,  'LPt ' ,0, 'NEW' , 'NC' ) 

CALL  CLS 
10  TYPE  500 

00  FORMAT  (////20X,'G  R  A  P  H  I  N  G '/// 1  OX ,' PAPER  READY  Y/N  ? ' / ) 
READ  (5,505)  PAP 
05  FORMAT  (Al) 

IF  ( P A P ♦ N E .  ' Y '  )  GO  TO  510 
TYPE  520 

20  FORMAT  (10X,'PEN  READY  Y/N  ?'/) 

READ  (5,525)  PEN 
25  FORMAT  (Al) 

IF  ( PEN . NE . ' Y'  )  GO  TO  510 
30  TYPE  540 

40  FORMAT  ( 1  OX ,' PLOTTER  ACTIVE  Y/N  ?'/) 

READ  (5,545)  PLT 
45  FORMAT  (Al) 

IF  (PLT . N E » ' Y ' )  GOTO  530 

DRAW,  PLOT,  ETC 


o  o  n  o  on 
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C  ACTIVATE  PLOTTER 

WRITE  ( 3  . 650 ) 

PLOTTER  ACTIVE.  HOME(BOT,L/H  CNR  )  .  POSITION  2.CM  UP  AND 
ALONG. ORIGIN  AT  THAT  POINT 
650  FORMATS  H  A100.200  0  ') 


GO  TO  888 

DRAW  BORDER  OF  GRAPH 
590  FORMAT! 'R' . 15. . 15. '  '> 

595  FORMAT  < ' M20  ') 

C  595  =  TICK  FOR  X-AXES  AND  Y-AXES 

WRITE  (3.565) 

DO  600  1=1 >15 
WRITE  (3.595) 

C  WRITE  (3.565) 

WRITE  (3.605) 

605  FORMAT ( 'R150.0  '  ) 

600  CONTINUE 

C 

DO  610  1=1.6 
WRITE  (3.595) 

C  WRITE  ( 3.565) 

WRITE  (3.615) 

615  FORMAT ( ' RO . 200  ') 

610  CONTINUE 

C 

DO  620  1=1.15 
WRITE  (3.595) 

C  WRITE  (3.565) 

WRITE  (3.625) 

625  FORMAT ( 'R- 150.0  ') 

620  CONTINUE 

C 

DO  630  1=1.6 
WRITE  (3.595) 

C  WRITE  (3.565) 

WRITE  (3.635) 

635  FORMAT! 'RO. -200  ') 

630  CONTINUE 

WRITE  (3.575) 

C 

C  END  OF  BORDER 

n 

C  WRITE  LABELLING 

C 

L  A  B  =  0 
I XX=-60 
I Y  Y  =  -  5  0 

WRITE  (3.560)  IXX.IYY 
C  POSITION  ABSOLUTE 

DO  645  1=1.15 
WRITE  (3.655)  LAB 
655  FORMAT ( ' S12  ' » 1 2 . ' _ ' ) 

C  HEIGHT=2»5mm.R0TATI0N=N0RMAL 

L A  B  =  L  A  B  +  1 
WRITE  (3.705) 

C  MOVE  TO  NEXT  LOCATION 

645  CONTINUE 

WRITE  (3.655)  LAB 


CJ  CJ  CJ  CJ  O  L)  U  CJ  U  (J  u  u 
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705  FORMAT ( '  R 1 0  5 . 0  ') 

710  FORMAT (  '  R  ~  5  0  ?  0  ' ) 

C 

l_  A  B  =  0 
I XX=-80 
I Y  Y  =  -  2  5 

WRITE  (3.560)  IXX»IYY 
DO  660  1  =  1 . 7 
WRITE  (3.655)  LAB 
L AB=LAB+ 1 
WRITE  (3.615) 

WRITE  (3.710) 

C  STEP  VERTICAL 

660  CONTINUE 

C 

ANNOTATION 

I XX=650 
IYY=-120 

WRITE  (3.560)  IXX.IYY 
WRITE  (3.665) 

665  FORMAT ('SI  2  WAVELENGTH  ( u  M >  _  '  ) 

H  E I G  H  T  =  3  m  m .  R  0  T  A  T 1 0  N  =  N  0  R  M  A  L 
1XX=-100 
I YY=200 

WRITE  (3.560)  IXX.IYY 
WRITE  (3.670) 

670  FORMATS S42  EXTINCTION  COEFFICIENT  ( SQH/GM  )  _ '  ) 

HEIGHT=3 mm. ROTATION =270” 

IXX  =  0 
I Y  Y  =  1  2  3  0 

WRITE  (3.360)  IXX.IYY 
WRITE  (3.675)  TH.RH.CONC 

X 

673  FORMAT ( 'St 2  TEMP : ~ ' . F6 . 2 . ' C  REL ♦ HUMIDITY : -  '  .  F6 . 2 .  '  %  CONCE 

*NTRATION:-  '  »  F  6 . 4  »  '  GM/CUB»METRE_'  ) 

I X  X  =  4  5  0 
I Y  Y  =  1  3  5  0 

WRITE  (3.560)  IXX.IYY 
WRITE  (3.680)  DTI . DT2.DT3.ASMKT » I  RUN 
680  FORMAT ('SI  2  DATE:-  ' . A2 . ' / ' . A3  .  '  /  '  A2 ?  '  TIMEt-'.A4.'  RUN 

*  *  '  .  1 2  .  '  _  '  ) 

I XX=200 
I YY=1290 

WRITE  (3.560)  IXX.IYY 
WRITE  (3.685)  ( RING ( I ) . 1=1 . 10 ) 

685  FORMAT  <  'S12  SMOKE  TYPEJ-  '.10A4.'_') 

I X  X  =  .1 5  0 
IYY=1450 

WRITE  (3.560)  IXX.IYY 
WRITE  (3.690) 

690  FORMAT  ('S13  MRL..  SMOKE  CHAMBER  MEASUREMENTS..') 

END  OF  ANNOTATION 
560  FORMAT (  'A'  .  15.  '  .  '  .  15.  '  ') 

PLOTTING  ACTUAL  VALUES  FROM  DATA 


ou  uo 
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720  FORMAT  (  '  M  2 1  '  ) 

565  FORMAT (  '  D  ') 

C  PEN  DOWN 

575  FORMAT ('U  ') 

C  PEN  UP 

WRITE  ( 3  . 575  > 

C 

888  IXX=INT( ( ALAMDA(3>*150)+0.5> 
I Y Y  =  I  NT ( (ACALC(3>*200)+0,5> 
CALL  LOCLIN 
DO  570  1=3-100 

IXX=INT<  (  ALAMDA(I>*150)+0.5> 
IYY=INT( < ACALC( I ) *200) +0,5) 
CALL  LOCLIN 
570  CONTINUE 

WRITE  (3.575) 


IXX=INT< (ALAMD(3)#150>+0.5) 

I Y Y  =  I  NT ( (ACALC( 106 ) *200 ) +0 ♦ 5 > 

CALL  LOCLIN 

DO  580  1=106.202 

IXX=INT(  (  AL. AMD  (  1-103  )*  150)  +0*5) 

IYY--INT(  (  ACALC  (  I )  *-2  0  0  )  +  0 ,5) 

CALL  LOCLIN 
580  CONTINUE 

WRITE  (3.575) 

C 

IXX  =  INT<  (ALAM(3>*150)+0,5) 

I YY= I  NT ( ( ACALC  <  209  >  #200 ) +0 . 5 ) 
CALL  LOCLIN 
DO  585  1=209.305 
[ XX= I  NT ( ( AL AM ( I -206 ) *1 50 ) +0 ♦ 5 ) 

I Y  Y  = I N  T (  (ACALC ( I) #200) +0,5) 

CALL  LOCLIN 
585  CONTINUE 

WRITE  (3.575) 

C  PEN  UP  -  AT  END  OF  PLOTTING 

r 

■  W  R I T  E  (3.640) 

640  FORMAT ("H  0  ') 

C  INACTIVATE  PLOTTER 

700  FORMAT! ' A' . 16.  '  .  '  .  16.  '  '  ) 

CALL  CLOSE  (3) 


TYPE  550 

550  FORMAT  (///lOX.'END  OF  GRAPHING 10X PRESS  RETURN  KEY'/) 
READ  (5.555)  DUD 
555  FORMAT  (Al) 

RETURN 

END 

C 

C 

C  LOCATES  START  OF  LINE  AND  PUTS  PEN  DOWN 

C 

SUBROUTINE  LOCLIN 


C 
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*ALAMD( 103) * ALAM( 105) ,RING( 10) ,DT1 *  DT2  * DT3  *  AHSKT  * INSKR  *  ASMKT 
*  * TH  *RH,C0NC* IBAKGD* IRUN  , NRUN ( 50 ) *  I X  X  *  I Y  Y  * I  BACK (312) 

IF  < IYY.GT. 1200. OR. IYY.LE.O)  GO  TO  335 
WRITE  ( 3  * 300 )  IXX*IYY 
WRITE  (3*305) 

C  PEN  DOWN 

330  RETURN 

335  WRITE  (3*310) 

C  PEN  UP 

GO  TO  330 

300  FORMAT ( 'A' *16* ' * ' ,16*  '  ') 

305  FORMATED  ') 

310  FORMAT (  '  U  '  ) 

END 


r 


ISMKR 


ERL- 0386 -MA 


46 


APPENDIX  V 


MRLV21.FOR 


C 


c 


c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

r* 

L- 


c. 


c 


c 


c 


700 

795 

705 

710 


711 

C 


715 

716 


C 


720 


721 


FILENAME:-  MRLV2 1  *  FOR 

old  filename:-  MRL 1  *  FOR 

LINK  OB JS  OF  MRLV20.F0R»MRLV21.F0R»MRLV22.F0R»READ43.MAC 

MRL  SMOKE  CHAMBER  PROGRAMME 
ACCEPTS  GOLAY  CELL  INPUT 
STORES  S  ANALIZES  DATA 

MODIFIED  TO  USE  312  BACKGROUND  READINGS 

RE-RUN  OF  EXISTING  DATA  STORED  ON  DISC 
SUBROUTINE  RERUN 

COMMON  I  TEMP (2200) » INOSMK < 312 )» ISMK < 312 ) ?ACALC(312) r ALAMDAC  103)  , 
*ALAMD< 103) »  A L A M ( 105) » RING ( 10) » DTI » DT2 » DT3 » ANSKT » INSKR  t ASMKT » ISMKR 
*  t  TH  ?  RH  *  CONC  ?  I BAKGD  r  I  RUN  ?  NRIJN  (  50  )  ?  IXX  ?  I YY  t  I  BACK  (312) 


INTEGER *2  DINK ( 7  > ? DISK ( 7 ) ? DICK (7 ) 


DATA  D INK/ ' DY ' » ' 1 J ' » ' IN  '  »  '  SK  '  »  '  00 ' >  '  ♦  D  '  f 
DATA  D I S  K  /  '  D  Y  '  ?  '  1  :  '  ?  '  I S  '  ?  '  M  K  '  ?  '  0  0  '  ?  '  ♦  D  '  ? 
DATA  D I CK/ ' DY  '  »  ' 1 :  ' t  ' I C ' ? ' LK  '  >  ' 00  '  t  '  * D ' ? 


AT  '  / 
AT'  / 
AT  '/ 


TYPE  700 

FORMAT  (10X)'W  RITE  DOWN  EX 
#M  B  E  R'/10Xf'N0W  PRESS  RETURN  KEY') 
READ  <  5  ?  7  9  5  >  DUD 
FORMAT  ( A  1  ) 

CALL  CLS 


STING  RUN 


TYPE  710 

FORMAT  (1  OX?  "RE-RUN 
*R  U  N  N  U  M  B  E  R 
READ  (5>7li)  I RUN 
FORMAT  (12) 


OF  PREVIOUS  DISC-STORED  DATA ' //5X »' GIVE 
YOU  WANT  TO  SEE') 


N  U 


D I N K<5)-NRU N ( I R  U N  ) 

CALL  ASSIGN  < 2 » DINK » 1 4 » ' OLD  '  > 

READ  ( 2  »  7 1 5  )  ( INOSMK ( I ) j 1  =  1 » 312 ) 

FORMAT  <63<5< I10?2X) »/>  ) 

READ  ( 2  »  7 1 6  )  DT 1 ?  DT2 ?  DT3  ? ANSKT  j INSKR 
FORMAT  (  A2  ?  A3  ?  A2  ?  5X  ?  A4  ?  5X?.I2  > 

CALL  CLOSE  (2) 

DISK(5)=NRUN( I  RUN) 

CALL  ASSIGN  (  2. ?  D  I  3K  .  1  4  ?  ' OLD  '  ) 

READ  ( 2  j  720  )  (ISMK(I) ?  I  =  1  ?  3 1 2 ) 

FORMAT  ( 6 3 ( 5 (  1 1 0  ?  2 X ) 7 / )  ) 

READ  ( 2  t 721 )  DTI ?  D  T  2  ?  D  T  3  j  A  S  M  K  T ? ISMKR? ( R I N  G ( I ) j 1  =  1? 10) 
FORMAT  (  A  2  ?  A  3  ?  A  2  ?  5  X  *  A  4  ?  5  X  ?  1 2  » 1 0  A  4  ) 

CALL  CLOSE  (2) 


DICK(5)=NRUN( I  RUN) 

CALL  ASSIGN  C 2 ? D I CK ? 1 4 ?  '  OLD ' ) 

READ  ( 2  j  7 2 5  )  ( ACALC ( I  )  ? I - 1 »  3 1 2 ) 

FORMAT  ( 6 3 ( 5 ( F 6  *  3  ?  2 X  )  * /  )  ) 

RE!AD  (2?  726)  DT 1 ? DT2 ? DT 3 ? ASMKT ? CONC ? TH » RH ? I BAKGB 


n  o 
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726  FORMAT  < A2 » A3 » A2 , 5X » A4 » SX , F6 . 4 , 5X , FA ♦ 2 » 5X » F6 . 2 , 5X » 1 6 ) 
CALL  CLOSE  <2) 

C 


TYPE  730, IRUN» DTI ,DT2,DT3, ( RING ( I ) , 1  =  1 , 10) ,CONC 
730  FORMAT  ( / / 5  X , ' R  U  N  # ' ,  1 2 , 5  X  , 'BATE!-' ,  A  2  ,  A  3  ,  A  2  /  5  X  , "SMOKE  TYPE;-', 
♦10A4/5X, 'CONCENTRATION  ' , F 6 . 4 ,  '  GMS/CUB * METRE  ' //I  OX ,  '  I  R  THIS  THE 
♦CORRECT  RUN  ?  Y/N') 

READ  <5,735)  COR 


735 


740 

745 


FORMAT  ( A 1 ) 

IF  ( COR . EQ . ' Y ' )  GO  TO  740 
GO  TO  705 
TYPE  745 

FORMAT  ( / / 2 0 X ,  ' R  E  M  E  M  B  E  R  !  !' /5X ,' RE-INSERT  WRITTEN  DOWN  R 

♦UN  NUMBER  BEFORE  CLOSING  DOWN  COMPUTER ,  ' /5X ,  ' OR  CARRYING  ON  WITH  A 
♦NOTHER  LEGITIMATE  RUN , ' //5X , ' OR  THE  DISC  DATA  FILES  MAY  BECOME  OVE 
♦  R  —  WRITTEN  !  !'///5X,' PRESS  RETURN  TO  GET  TO  M  E  N  U ' ) 

READ  <5,795)  DUD 

RETURN 

END 


C 


C 

C 


C 

C 

C 

c 


♦ 

♦ 

♦ 


SUBROUTINE  ERROR < V 1 , 02 , I  ERR ) 


DIMENSION  UV2  ( 17) , AV2( 17) 

D h T A  U V 2 / *0155, *024, *0335,  *044,  *054, *064,  * 
.24, ,332, .44, .54, .64, .738, *825, .92/ 

DATA  A 0 2 / . 0 0 3 9 ,  .0027,  *00223,  .00194,  .00179, 
.00153, .0013, *00118,  ,00111 , .00105, *000102, 
.00096, .00094/ 


074, .084, .014, 

*00167, *0016, 
.00099, .00097, 


I E  R  R  =  0 

IF  (VI. LT, 0.025)  GO  TO  1 


IF(Vl.LT.O.l)  GO  TO  10 


♦♦♦♦♦♦ 


VALUES  OF  VI  BETWEEN  0,1  AND  1,0 


DO  29  .1-1,9 

A - 0 . 1+0. 1 ♦FLOAT <1-17 
B  =  A  +0 , 1 

I F ( V 1 .GE.A.AND.Vl ,LT  *B)  GO  TO  21 
GO  TO  29 

21  I F ( V 2 . L T . A V 2 ( J ) . 0 R . V 2  *  G T  *  U V 2  <  J )  )  GO  TO  1 

GO  TO  99 
29  CONTINUE 
GO  TO  99 

**#*♦♦♦•#♦♦♦♦♦♦♦♦***♦♦♦♦♦♦*•♦♦*♦*♦*♦♦♦♦*♦*  **♦♦*  ♦*#**♦**♦♦#♦♦ 


C  VALUES  OF  VI  BETWEEN  0.025  AND  0*1 

C 

10  IF < VI. GE. 0,025. AND.V1.LT. 0,03)  GO  TO  11 
DO  19  1  =  1 ,7 

A  =  0. 03  +  0, 01 ♦FLOAT  < I - 1  ) 

B  =  A  +  0 , 01 

IF<V1 . GE.A.AND.Vl. LT.B)  GO  TO  12 
GO  TO  19 

I F ( V 2 , L T , A V 2  < I  + 1 ) , 0 R , V 2 , G T . U V 2 (  I  +  1  )  )  GO  TO  1 
GO  TO  99 


12 
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1?  CONTINUE 
GO  TO  99 

**********************  jit*******#*)!!*****)!;*************;******* 
VALUES  OF  01  BETWEEN  0.025  AND  0.03 
UR  I TE ( 3  f  # ) 0 1 >02 

11  I F ( 02 . LT . A02 ( 1 ) . OR ♦ 02 . GT  . U02 ( 1 )  )  GO  TO  1 
GO  TO  99 
C 

C  *********************************************************** 

c 

C  DATA  UNACCEPTABLE  DUE  TO  UNACCEPTABLE  ERRORS 
C 

1  I ERR= 1 

99  CONTINUE 

RETURN 
END 
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APPENDIX  VI 
MRLV22 .FOR 


C 

C 

C 

C 

C 

C 

C 

C 

C 

C 


C 

c 

c 

c 

c 


c 


c 

c 

c 

c 

c 


c 


soo 


805 


310 

830 

840 

345 

850 

820 

C 

C 

C 


FILENAME:-  MRLV22.F0R 

OLD  FILENAME  J-  MRL2.F0R 

LINK  OB JS  OF  MRLV20.F0R»MRLV21 . FOR » MRLV22 . FOR » READ43.MAC 
MODIFIED  TO  USE  312  BACKGROUND  READINGS 


SUBROUTINE  CL.S 
CLEARS  0  D  U  SCREEN 
TYPE  5 

FORMAT  <//////////////////////////////) 

RETURN 

END 


INITIALLIZE  IRUN  COUNTER  FOR  NEW  DISKETTE 
SUBROUTINE  WEN 


COMMON  ITEM P (2200) t IN0SMK(312) > 1 SMK ( 3 1 2 ) *ACALC(312) t ALAND A ( 103 ) , 
♦AL AMD  (.103)  t  ALAH  ( 105)  *  R  I N G  (  1 0  )  t  DTI ,  DT2  t  DT3  >  ANSKT  t  INSKR  t  ASMKT  t  ISMKR 
♦»TH>RH>CONC» IBAKGDj I  RUN t NRUN ( 50 ) ? I XX  *  I Y  Y  « IBACK(312) 


I RUN= 1 
RETURN 
END 


INCREMENT  TRIAL  RUN  NUMBER 
SUBROUTINE  NUR 

COMMON  I  TEMP ( 2200 ) » INOSMK( 312) t ISMK ( 3 1 2 ) » AC ALC ( 3 1 2 ) t ALAMDA ( 103) t 

♦  ALAND  (  103)  ,  ALAM(  105)  »  R I N  0  ( 10)  >  DTI  «DT2»DT3.»  ANSKT  *  INSKR  t  ASMKT  »  ISMKR 
*t TH » RH»CONC» IBAKGD> IRUN » NRUN (50 ) . I XX  t I Y Y  ? I  BACK (312) 

CALL  CLS 
TYPE  800 

FORMAT  ( 1  OX  f  ' CHANGE  RUN  NUMBER  ' //I OX j  ' DO  YOU  WANT  TO  INCR 

♦  EMENT  RUN  NUMBER  t  ' / '  OR  CHANGE  THE  RUN  NUMBER  ?  TYPE  I/C') 

READ  ( 5 i 805 )  CH 
FORMAT  (Al) 

IF  ( CH ♦ EQ ♦ ' C  '  )  GO  TO  S10 
IRUN-IRUN+1 
GO  TO  820 
TYPE  330 

FORMAT  ( 1 0 X  f  'GIVE  REQUIRED  RUN  NUMBER'  ) 

READ  ( 5  r S40  )  IRUN 
FORMAT  (12) 

TYPE  345 r  IRUN 

FORMAT  ( 2  0  X  ?  'RUN  *  , I2//10X , ' PRESS  RETURN  KEY') 

READ  ( 5  >  350 )  DUD 
FORMAT  (Al) 

RETURN 

END 
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C 

C 


C 

C 

C 

C 


C 

C 

C 

C 

C 

C  200 

c 

270 


275 


280 

214 

C 

211 


210 

220 

■~>ncr 

uj 

222 

212 

230 

231 
C 

C 

C 


221 

21? 


NO  SMOKE  &  BACKGROUND  SUBROUTINES 
SUBROUTINE  NSK(.IBAK) 

COMMON  I  TEMP (2200 ) . I NOSMK  (312). I SMK (312) *  ACALC (312)?  ALAMDA (103) . 
*ALAMD( 103) » ALAM( 105) >RING( 10) . DTI .DT2.BT3. ANSKT. TNSKR. ASMKT . ISHKR 
#»TH» RH»C0NC»IBAKGD»IRUN»NRUN(50) , I XX, I YY . I B ACK ( 3 1 2 ) 


1NTEGER*2  D I SK ( 7 ) . D I S ( 7 ) 


DATA  DISK/'DY'  »  ' 1 ! S ' IN' ,  'SK ' » '00' .  '  .D'  .  'AT'/ 
DATA  DIS/'DY ' . ' 1 l ' t  ' IB' , ' AK' . '00 ' , ' ,D' . ' AT'/ 


CLEAR  NO  SMOKE  DATA  FIELD 
DO  200  I  =  1.312 
INOSMK(I)  =  0 
CONTINUE 

NO  SMOKE  DATA  FIELD  CLEARED 
TYPE  270 

FORMAT  ( // .  3X »  'USE  EXISTING  DATA  (NOT  MORE  THAN 
*  ONE  HOUR  OLD)  Y/N') 

READ  ( 5  >  275  )  Y 
FORMAT  (Al) 

IF  <  Y ♦ NE ♦ ' Y  '  )  GO  TO  280 
IF  ( I B  A  K  .  E  Q  .  0  )  GO  TO  213 
GO  TO  223 
CALL  CLS 

IF  (IBAK.EG.O)  GO  TO  210 
TYPE  211 

X 

FORMAT  (///YIOX.'B  A  C  K  6  R  0  U  N  D  R  E  A  D  I  N  6  S'//l 
*0X.'B  LANK  OFF  N  E  R  N  S  T  GLOWER'/) 

GO  TO  212 
TYPE  220 

FORMAT  (  /  /  /  /  2  0  X  »  '  N  0  SMOKE  M  E  A  S  U  R  E  M  E  N  TV) 

TYPE  225 

FORMAT  (20X»'GIVE  TIME  EG  0900') 

ACCEPT  222  .ANSKT 
FORMAT  ( A  4 ) 

TYPE  230 

FORMAT  ( / 1  OX  > 'GIVE  A  NUMBER  BETWEEN  1  i  7  FOR  REVS  REQUIRED'/) 
READ  (5.231)  INSKR 
FORMAT  (12) 

IF  ( INSKR.LT. 1 .OR. INSKR.GT.7)  GO  TO  214 

CALL  COLDAT ( ITEMP. INSKR ) 

IF  (IBAK.EG.O)  GO  TO  21S 
AVERAGE  BACKGROUND  READINGS 
N=  I NSKR 

DO  21?  1-1.311 
IB=0 

DO  221  J  =  1»N 

IB  =  IB  +  ITEMF‘(311*(  J-l  )  +  I  ) 

CONTINUE 
I  BACK (  I  )  =  I R/N 
CONTINUE 
GO  TO  223 
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C 

C  AVERAGE  NOSMOKE  DATA  AND  STORE  IN  NSMK  ARRAY 

218  N  =  INSKR 

DO  240  I  =•  1,311 
IB  =  0 

DO  250  J  =  1,N 

IB  *  I B  + 1  TEMP  ( 31 1  * ( J- 1 ) +1 ) 

250  CONTINUE 

INOSMK (  I  )  =  IB/N-IBACK<  I  ) 

240  CONTINUE 
C 
C 
C 
C 

C  STORE  ON  DISC  ?  ***USE  CURRENT  RUN  NUMBER 

C  AS  PART  OF  FILE  NAME 

C 

223  IF  (IBAK.EQ.O)  GO  TO  213 
D I S  <  5 )  =  NRUN(IRUN) 

CALL  ASSIGN  ( 2 > D I S , 1 4 , ' NEW ' > 

WRITE  (2,243)  ( I  BACK < I ) , I  =  1 , 3 1 2 ) 

CALL  CLOSE  (2) 

IF  ( Y.EQ. ' Y'  )  GO  TO  266 
TYPE  216 
GO  TO  217 

C  X 

216  FORMAT  (5X,'R  EMOVE  SOURCE  BLAN  K'/10X,'T 

♦  HEN  PRESS  ENTER  KEY') 

213  D I S  K  (  5  )  =  NRUN(IRIJN) 

CALL  ASSIGN  (2, DISK, 14, 'NEW' ) 

WRITE  (2,243)  ( INOSMK (I), 1-1,3 12) 

243  FORMAT  (63(5(110, 2X),/)) 

WRITE  (2,244)  DTI , DT2 , DT3 , ANSKT , INSKR 

244  FORMAT  (A2,A3,A2,5X,A4,5X, 12) 

CALL  CLOSE  (2) 

IF  ( Y  ,  EG  .  ' Y  '  )  GO  TO  266 
C 

C  DATE, NOSMOKE  START  TIME  AND  NUMBER  OF  REVS  -ON  END  OF  DATA 

C 

217  TYPE  260 

C  X 

260  FORMAT  (//10X,'END  OF  MEASUREMENT ', 1  OX , ' PRESS  RETURN  KEY'/) 

READ  (5,265)  DUD 

265  FORMAT  (Al) 

266  RETURN 
END 

C 

C 

C 

C  SMOKE  SUBROUTINE 

C 

SUBROUTINE  SK 
C 

COMMON  I  TEMP (2200) , INOSMK (312) , ISMK(312) ,ACALC(312) , A L A M D A ( 103) , 

♦  ALAMO (103) ,ALAM( 105) , R I  NO ( 10 ) , DT 1 , DT2 , DT3 , ANSKT , I NSKR , ASMKT , I SMKR 
*,TH,RH,CONC» IBAKGD, I  RUN , NRUN ( 50 ) , I XX, IYY, IBACK(312) 

C 

INTEGER *2  DISK(7) 

C 

C 
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DATA  DISK/'DY' * ' 1 J ' * ' I S ' * ' MK ' * ' 00 ' , ' . D  '  * ' AT ' / 


DATA  RING/10*'  '/ 


CLEAR  DATA  FIELD 
DO  300  I  =  1 f  3 1 2 
I SMK  < I )  =  0 
300  CONTINUE 

SMOKE  DATA  FIELD  CLEARED 
CALL  CLS 
310  TYPE  320 

320  FORMAT  <////20X*'S  MOKE  MEASUREMENT'/) 

TYPE  325 

325  FORMAT  <20X*'GIVE  TIME  EG  0930') 

ACCEPT  322-ASMKT 
322  FORMAT  (A4) 

TYPE  360 

360  FORMAT  ( / 1  OX , ' G I VE  SMOKE  TYPE/ DESCRIPTION'/) 

ACCEPT  370* (RING(I) >1=1*10) 

370  FORMAT  (10A4) 

TYPE  330 

380  FORMAT  (/10X»'GIVE  CHAMBER  TEMPERATURE  IN  DEGREES  CENT.'/) 
ACCEPT  *  *  TH 
TYPE  375 

375  FORMAT  (/10X*'GIVE  RELATIVE  HUMIDITY  IN  PERCENT'/) 

ACCEPT  *  *  RH 
TYPE  230 

230  FORMAT  </9X*'GIVE  A  NUMBER  BETWEEN  1  8  7  FOR  REA'S  REQUIRED' 
READ  (5*332)  ISMKR 
382  FORMAT  (12) 

IF  (ISMKR. LT.l, OR. ISMKR. GT.7)  GO  TO  310 
CALL  COLDAT( ITEMP* ISMKR) 

AVERAGE  SMOKE  DATA  AND  STORE  IN  SMK  ARRAY 
N  =  ISMKR 
DO  330  I  =  1*311 
IB  =  0 

DO  340  J  =  1*N 

IB  =  I B  + 1 T E M P ( 3 1 1  * ( J  - 1 )  + 1  ) 

340  CONTINUE 

ISMK(I)  =  IB/N-IBACK ( I ) 

330  CONTINUE 
C 
C 
C 

C  STORE  ON  DISC  ?  ***USE  CURRENT  RUN  NUMBER 

C  AS  PART  OF  FILE  NAME 

C 

DISK (5)  *  NRUN(IRUN) 

CALL  ASSIGN  (2* DISK *14* ' NEW ' ) 

WRITE  (2*342)  ( I SMK ( I ) *  I  =  1  *  3 1 2 ) 

342  FORMAT  ( 63 ( 5 ( 1 1 0  * 2X ) * / )  ) 

WRITE  (2*344)  DTI  * DT2 * DT3 * ASMKT * ISMKR * ( RING ( I ) *  1  =  1  *  10 ) 

344  FORMAT  ( A2 *  A3  * A2 * 5X * A4 * 5X *  12  *  1 0A4 ) 

CALL  CLOSE  (2) 

C 
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C 

C  DATE  *  SMOKE  START  TIME*  SMOKE  REVS  ON  END  OF  DATA 

C 

C 

C  X 

TYPE  350 

350  FORMAT  (//10X*'END  OF  SMOKE  MEASUREMENT' *  1 0  X  > 'PRESS  RETURN  KEY'/) 
READ  (5*355)  DUD 
355  FORMAT  (Al) 

RETURN 

END 

C 

C 

C 

C  CALCULATE  SMOKE  S  NOMSOKE  READINGS*  FOR  RELATIVE  VALUES 

C 

SUBROUTINE  CLC 
C 

COMMON  I  TEMP ( 2200 ) , INOSMK< 312) , ISMK(312) * ACALCC312) *  ALAMDA (103) * 
*ALAMD(  103)  *ALAM(  105)  *RING(10)  ,  DT.I  *  DT2  *  DT3  *  ANSKT  *  I NSKR  *  ASMKT  *  I SMKR 
#»TH»'RH*CONC»  IBAKBD*  I  RUN  *  NRUN  (  50  )  *  I  XX  *  I Y  Y  *  I  BACK  (  3 1 2  ) 


INTEGERS  D I  S K  (  7  ) 


DATA  DISK/ ' DY ' *  '  1  ;  '  >  '  IC  '  *  'LK ' * ' 00  '  *  '  . D ' >  '  AT  '  / 


CLEAR  CALC  DATA  STORE 
DO  400  1=1*312 
ACALC ( I ) =0 
400  CONTINUE 

CALC  FIELD  CLEARED 
CALL  CLS 
TYPE  410 

410  FORMAT  (////20X*'C  ALCULATE  DAT  A'/) 

IB  =  0 

IB=IB+INOSMK( 102) 

IB=IB+IN0SMK(205) 

IB=IB+IN0SMK(308) 

IB=IB+ISMK( 102) 

IB=IB+ISMK(205) 

IB=IB+ISMK(308) 

I B= I B/S 
IB  =  0 

TAKEOUT  I B — 0  *  ifc  $  $  *  * 

IBAKGD=IB 

C  BACKGROUND  AVERAGED 

V0L=4*5*2, 63*2*74 

C  VOLUME  OF  CHAMBER 

A  L  E  N  =  5  *  0 

C  LENGTH  OF  DETECTOR  PATH 

415  TYPE  420 

420  FORMAT  (//IOXj'GIVE  MASS  OF  SMOKE  USED  IN  GRAMS'/) 

READ  (5**)  AMASS 
CONC=AMASS/VOL 
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DO  430  1=1,312 

IF  ( I NOSMK  < I ) . EQ . 0 )  GO  TO  430 
IF  (ISMK(I) .  E  Q  .  0  )  GO  TO  430 
G0LA=FL0AT ( I SMK ( I ) ) /FLOAT ( I NOSHK (  I  )  ) 

IF  ( GOLA . LE  *  0 ♦ 0  )  GO  TO  430 
ACALC( I )=(-AL0G(GQLA) ) / ( ALEN*CQNC ) 

C  %%%%&%%%){%%%%%%%%%%  %%%%%%%%%  %%% 

VI  =  (  FLOAT*:  I  NOSMK  (  I  )  )  )/2000 
V2-(FL0AT  < I S  M  K  ( I ) ) )/2000 
CALL  ERROP  (V1,V2,IERR> 

IF  (IERR.NE.O)  ACALC ( I ) =0 . 0 

430  CONTINUE 
C 
C 
C 

C  STORE  ON  DISC  ?  ***USE  CURRENT  RUN  NUMBER 

C  AS  PART  OF  FILE  NAME 

C 

D I S  K ( 5 )  =  NRUN(IRUN) 

CALL  ASSIGN  ( 2 , DISK  ,  1 4 , ' NEW  '  > 

WRITE  (2 >442)  ( ACALC (I), 1=1,312) 

442  FORMAT  < 63 ( 5 ( F6 . 3 , 2X ) , /  )  ) 

WRITE  ( 2  j  44 4 )  DT 1 , DT2 , DT3 , ASMKT  , CONC , TH , RH , I BAKGD 
444  FORMAT  ( A2 , A3 f A2 j 3X t A4 r 5X » F6 . 4 » 5X t FA . 2 t 5X » F6  ♦  2 > 5X » 1 6  ) 
CALL  CLOSE  (2) 


DATE, SMOKE  START  TIME,  CONCENTRATION , TEMPERATURE 
RELATIVE  HUMIDITY, BACKGROUND  READING--  ON  END  OF  DATA 


TYPE  440 

440  FORMAT  <//10X,'END  OF  CALCULATIONS' ,10X, 'PRESS  RETURN  KEY'/) 
READ  (5,445)  DUD 
445  FORMAT  (Al) 

RETURN 

END 


TABULATE  RESULTS  SUB 
SUBROUTINE  TABUL 


COMMON  I  TEMP < 2200 ) , I NOSMK ( 312) , I SMK (312) , AC ALC ( 3 1 2 ) , AL AMD A ( 1 03 ) , 
*ALAMD(103) ,ALAM( 105)  , RING ( 10) , DTI , DT2 , DT3 , ANSKT  ,  INSKR , ASMKT , ISMKR 
*,TH,RH,CONC, IBAKGD, IRUN,NRUN(50) , I XX, IYY, I  BACK (312) 

C 

C 

C 


CALL  CLS 

TYPE  600 

AOO 

FORMAT  ( / / / / 2 0 X ,  '  T  A 

B  U  L. 

ATE  RESULT  S'//) 

615 

TYPE  605 

605 

FORMAT  ( 10X, 'PRINTER 
READ  (3,610)  PR 

READY 

Y/N  ?') 

610 

FORMAT  (Al) 

no  on  o  o  no 
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IF  (PR.NE.  '  Y  '  )  GO  TO  615 
C 
C 

TYPE  620 

C  PRINT  620 

620  FORMAT  (//////) 

TYPE  625 

C  PRINT  625 

625  FORMAT  C//7X?'M  RL.  SMOKE  CHAMBER  MEASUREM 

1 ============'//) 

C  PRINT  630? IRON, DTI ?DT2? DT3 

TYPE  630? IRUNfDTl ? D  T  2  t  D  T  3 

630  FORMAT  C10X?'RUN  * ' ? 1 2 ? 20X ? ' DATE J -  ' ? A2 ? ' / ' ? A3 ? ' / ' ? A2/ ) 

C  PRINT  635? ANSKT? INSKR 

TYPE  635?ANSKT  ? INSKR 

635  FORMAT  ( 3X ?' NOSMOKE : -  TIME  =  '  » A4 » 3X »  ' RE VS  =  ',I2> 

C  PRINT  637? ASMKTf ISMKR 

TYPE  637? ASMKT? ISMKR 

637  FORMAT  < 5X ?' SMOKE ! -  TIME  =  ' ? A4 ? 3X ? ' REVS  =  '?I2/) 

C  PRINT  640?TH?RH 

TYPE  640  ?  TH  ?  RH 

640  FORMAT  (  3X  ?  ?  CHAMBER  TEMPERATURE*.-  '?F6.2?'  C'?l  OX? 'RELATIVE  HUMID 


1ITYJ-  '  ?  F6 . 2  ? ' % ' ) 

X 

PRINT  645? (RING( I ) ? 1=1 ? 10) 

TYPE  645? ( R I NG ( I ) ? 1=1 ? 10) 

645  FORMAT  <3X? 'SMOKE  DESCRIPTION!-  '?10A4) 

PRINT  642  ? CONC 
TYPE  642? CONC 

642  FORMAT  < 1 5X» ' CONCENTRATION J -  '?F6.4?  '  GMS/CUB . METRE  '  ) 

PRINT  647  ? IBAKGD 
TYPE  647? IBAKGD 

647  FORMAT  ( 1 5X ?' BACKGROUND  READING J » 16/) 

PRINT  650 
TYPE  650 

X 

650  FORMAT  <3X?'FILTER  *'?6X?'EXT  COEFF ' ? 6X ? ' L AMD A  UM '? 6X ?' NOSMOKE  * 

1 ' ?6X?  'SMOKE  *'  ) 


J=0 

DO  660  I=l?100 
K=I+J 

C  PRINT  665?K?ACALC< I ) ?  A  L  A  M  D  A ( I ) ? INOSMKC I ) ? ISMK< I  ) 

TYPE  665  ?  K  ?  AC ALC  < I ) ?  A  L  A  M  D  A ( I ) ? INOSMKC  I  )  ?  ISMKC  I  ) 

665  FORMAT  ( 4X ? 1 5 ? ?X ? F6 . 3 ? SX ? F6 . 3 ? 8X ? 1 6 ? 8X ? I  6 ) 

660  CONTINUE 

C 

J  =  -3 

DO  670  I  =  1 04  ?  203 
K  =  I  +  J 

C  PRINT  665?K?ACALC(I)?ALAMD(I-103)?IN0SMK(I)?ISMK(I) 

TYPE  665?K?  ACAL.CC  I  )  ?  ALAMO  (  1-103  )  ?  INOSMKC  I  )  ?  ISMKC  I  ) 
670  CONTINUE 

C 

J  =  -6 

DO  675  1=207  ?306 
K  =  I  +  J 

C  PRINT  665?  K?  ACALCC  I  )?  ALAMC  1-206  )?  INOSMKC  I  )?  ISMKC  T. ) 


on  n  o  n  n  o  no  on 
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TYPE  665»K»ACALC(I) , ALAN (1-206) , IN0SMK( I)  ,  ISMK(  I  ) 
675  CONTINUE 

PRINT  620 
TYPE  620 


TYPE  690 

690  FORMAT  ( / 1  OX , 'TABLE  COMPLETED' ,  1 0  X  , 'PRESS  RETURN  KEY'/) 
READ  (5,695)  DUD 
695  FORMAT  ( A1 ) 

RETURN 

END 


FINISH  PROG  BY  STORING  ALL  RE LAMENT  DATA  ON  DISC 
SUBROUTINE  FIN 


COMMON  I  TEMP ( 2200 ) , IN0SMK(312) ,  ISMK(312) , ACALC ( 3 1 2 ) , ALAMDA < 103 ) t 
*ALAMD( 103) , ALAM( 103) , RING < 10) , DT 1 , DT2 , DT3 , ANSKT , INSKR , ASMKT , I SMKR 
#»TH»RH»CONC» IBAKGD, IRUN , NRUN ( 50 ) » I XX , I Y Y , I  BACK ( 312 ) 

C 

INTEGER*4  DISCO) 

DATA  DISC/ 'DY1  IRUN'  ,'  .DAT'/ 

C 

C 

IRUN=IRUN+1 

C 

CALL  ASSIGN  < 2 , D I  SC , 1 2 , ' NEW  '  ) 

WRITE  (2,888)  IRUN 
888  FORMAT  (12) 

CALL  CLOSE  (2) 

C 

RETURN 

END 
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Figure  2.  Nernst  Glower  characteristics 
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Figure  4.  Golay  cell  radiometer 
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Figure  5.  Smoke  IR  transmissometer  general  assembly 
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Figure  6.  Radiometer  rear  view 
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Figure  8.  Radiometer  electronics  block  diagram 
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Figure  9.  Radiometer  top  view 


Figure  10.  Golay  detector  (inside  cell  housing) 
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Figure  12.  Nichrome  on  glass  graticule  (read  by  HEDS1000  detector) 
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Figure  13.  HEDS1000  circuit  (comparator) 
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Figure  14.  HEDS1000  housing  and  circuit  board 
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Figure  17.  IR  synchronous  detector  circuit  (Card  1) 
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Figure  18.  Synchronous  detector  housing 


Figure  19.  Synchronous  detector  under  chassis  view 
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Figure  21.  Synchronous  detector  cardset 
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Figure  22.  Electronic  waveforms  timing 
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Figure  24. 


Golay  power  supply  and  preamp  box 


Figure  25.  Golay  power  supply  (Circuit  figure  23) 


Figure  26 


Figure  26 


Read  pulses  F10  circuit  +  Golay  preamp  (Circuit  figure  15) 
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Figure  27.  Data-logging/computer  block  diagram 
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Figure  28,  Empirical  graph  to  determine  aperature  material 
temperature  at  a  distance  from  a  hot  source 
(Nernst  Glower)  at  2000°C) 
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Figure  30.  Chopper  driver  circuit 
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Figure  32.  Nernst  Glower  AC  power  supply  and  chopper  driver 
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Figure  33.  Handheld  IR  beam  finding  detector 
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Figure  34.  Circular  variable  filter  simulator  circuit 
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Figure  35.  Subroutine  MRLGRP 
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Figure  36.  Subroutine  LOCLIN 
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Figure  37.  Subroutine  RERUN 
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Figure  38 


Figure  38.  Subroutine  CLS 
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Figure  39.  Subroutine  RAW(IFIL) 
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Figure  40.  Subroutine  MRLV20.FOR 
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Figure  41.  Subroutine  GRP 
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Figure  42.  Subroutine  RERUN 


ERL-0386-MA 
Figure  43 


Figure  43.  Subroutine  ERROP 
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N.B.  Initialise  IRUN  counter  for  new  diskette 


Figure  44.  Subroutine  WEN 
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Figure  45.  Subroutine  NUR 
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Figure  46,  Subroutine  NSK 
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Figure  47.  Subroutine  SK 
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Figure  48.  Subroutine  CLC 
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Figure  49.  Subroutine  TABUL 


Figure  50.  Subroutine  FIN 


ERL-0386-MA 
Figure  51 


Figure  51.  Macro  subroutine  COLDAT 
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Figure  52.  Subroutine  READ 
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Figure  53.  Subroutine  PRMT 


Figure  54.  Subroutine  ECD 
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Figure  55.  Subroutine  CNVT 
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Figure  56.  Subroutine  TRS 
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Figure  57.  Subroutine  TIM2 
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Figure  58.  Subroutine  TIM1 
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